ON THE CLASS 71SI OF RATIONAL SCHUR 
FUNCTIONS INTERTWINING SOLUTIONS OF 
LINEAR DIFFERENTIAL EQUATIONS 



DANIEL ALPAY, ANDREY MELNIKOV, AND VICTOR VINNIKOV 

Abstract. In this paper we extend and solve in the class of func- 
tions IZSX mentioned in the title, a number of problems originally 
set for the class IZS of rational functions contractive in the open 
right-half plane, and unitary on the imaginary line with respect 
to some preassigned self-adjoint matrix. The problems we con- 
sider include the Schur algorithm, the partial realization problem 
and the Nevanlinna-Pick interpolation problem. The arguments 
rely on the one-to-one correspondence between elements in a given 
subclass of TZST and elements in TZS. Another important tool in 
the arguments is a new result pertaining to the classical tangential 
Schur algorithm. 
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1. Introduction 

Functions iS'(A), which are C^^^- valued, analytic and contractive in the 
open right half plane C+, or equivalently, such that the kernel 

X + w* 

is positive in C+, play an important role in system theory, inverse 
scattering theory, network theory and related topics; see for instance 
[L] . |BC] . |DD] . [He], [X]. Here, positivity of the kernel means that for 
every n G N and every choice of points Wi, . . . ,Wn G C+ and vectors 
^i, . . . £ C^^^ the n X n Hermitian matrix 

is positive (that is, has all its eigenvalues greater or equal to 0). 

Far reaching generalizations of this class were introduced in [Ml IMVc] , 
in the study of 2D-linear systems (say, with respect to the variables 
(1^1)^2)); invariant with respect to the variable ti. To introduce the 
classes defined in these papers we first need a definition. 

Definition 1.1. Let ai, a2, 7 and'^^ he U^^^ -valued functions, contin- 
uous on an interval I = [a, 6]. Suppose moreover that and 02 take 
self-adjoint values, and that o\ is difjerentiable and invertible on I, and 
that the following relations hold: 

7(^2) + 7(^2)* = 7*(^2) + 7*(i2)* = -^cri{t2), t2 G I. 

dt2 

Then o"i, o"2, 7, 7* and the interval I are called vessel parameters. 

The class of functions SI corresponding to some vessel parameters cti, 
c"2; 7; 7* and I was introduced in [Ml IMVc] (see Definition 12.61 below) 
and consists of the functions 5'(A, ^2) of two variables A, t2 such that for 
every ^2 G I the function 5'(A, ^2) is meromorphic in C_|_ and the kernel 

Mh) - S{w,t2yai{t2)S{X,t2) 
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is positive for A and uj in the domain of analyticity of S{\,t2) in C+. 
For positive (11(^2), the positivity of the kernel imphes that S is analytic 
in C+; see [Do] . [X] . For general (invertible) cri(t2), the entries of S are 
of bounded type and S has only poles in C+; see [APRS] . It is also 
required that S'(A, ^2) is analytic at infinity for each ^2, with value Ip 
there, and that S'(A, ^2) niaps solutions of the input Linear Differential 
Equation (LDE) with the spectral parameter A 

d 

\(j2{t2)u{\ t2) - ai(t2)^«(A, ^2) + 7(^2)m(A, t2) = 0, 
to solutions of the output LDE 

d 

Xcr2*{t2)y{X,t2) - (yu{t2)^y{\t2) + 7*(^2)2/(A,t2) = 0. 

012 

It is proved in [Ml IMVc] that elements of SX are the transfer functions 
of ti-invariant conservative 2D systems; see Section [21 

The purpose of this paper is to study various questions in SI in the 
case of functions rational in A. A key result is the following theorem, 
which we prove in the sequel; see Section HI 

Theorem 14. II Le/^ us fix the parameters cti, a2, and 7, and the interval 
I. Then for every t^^ & \ there is a one-to-one correspondence between 
pairs {'j*,S) such that S G SX and 7* continuous in a neighborhood 
0/^2, and functions Y{X), meromorphic in C+, and with the following 
properties 

(1) r(oo) = Ip, 

(2) r(A)Vi(t°)r(A) < aiitl) for A G C+ where Y is analytic, and 

(3) Y{\)*ai{ti)Y{\) = for every A satisfying 3ftA = 0, and 
in a neighborhood of which Y is analytic. 

As mentioned above the cri(t2)-contractivity of Y implies that Y is 
of bounded type in C+, and thus the asserted non-tangential limits 
exist almost everywhere. But it is important to note that the theorem 
does not consider these limits, but only the points A in the domain of 
analyticity of Y. 



Definition 1.2. Functions Y with the properties in Theorem 4-1 will be 
called ai{t^) -inner functions, and their class will be denoted by S{t^). 
The subclass of rational functions of Sifl) will be denoted by T^Sif^) ■ 

For the sequel, it is important to notice that the general tangential 
Schur algorithm developped in [ADj can be applied to functions in 
Sifl), and in particular in IZS. 
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Definition 1.3. TZSX will denote the subclass of functions in SI 
which are rational in A for every t2 El. 

The paper consists of five sections besides the introduction, and we now 
describe their content. In Section [2] we review the main results from 
[M] and [MVc] on ti-invariant conservative 2Z)-systems, relevant to the 
present work. In particular the class SX mentioned above consists of 
the transfer functions of these systems. In Section [3] we present the 
reproducing kernel space approach to the tangential Schur algorithm 
for the class IZS, as developped in [ADj . We obtain in particular new 
formulas which allow us to find the main operator in a realization of an 
element of IZS after one iteration of the tangential Schur algorithm; 
see formulas fl3.15p . fl3.17p . (13.161) in Theorem 13. 5[ In Section |4] we 
develop the tangential Schur algorithm for a function S{X,t2) G IZSI. 
Applying directly the theory of the previous section to S{X,t2) leads 
to a new function which need not belong to IZSI. Instead, we ap- 
ply the Schur algorithm to the cri(t2)-inner function S{X,t2) for some 
preassigned ^2 ^ I? ^ind obtain a simpler (in terms of McMillan de- 
gree) (Ti(t2)-inner function So{X,t2). We use Theorem 14.11 to obtain 
an element in a class IZSI from 5'o(A,t2)- We call this procedure the 
tangential Schur algorithm for the class IZSI. We study in Section [5] 
the coefficients (called the Markov moments) Hi(t2) of the expansion 
of S{X,t2) around A = oo 



It turns out that the first Markov moment //o(^2) satisfies the Lyapunov 



7*(i2) - 7(^2) = (r2{t2)Ho{t2) - (ri{t2)Ho{t2)a^\t2)(r2{t2), G I, 



which means that given the functions o"i,o"2)7 and Hq on I, one can 
uniquely reconstruct 7*, and, as a result of Theorem 14.11 there will 
exists a unique function S{X, ^2) (^2 G I) with the given Markov param- 
eters. We solve the following problem: 

Problem 1.4 (Partial realization). Suppose that we are given the func- 
tions CTi, o"2 and the Markov moments Hi, i = 0, . . . ,n (n > 0) defined 
in 1, and satisfying there the (necessary) conditions (15. 5p 




i=0 




j=0 
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and (15.61) 

n— 1 

/or z = 0,...,n — 1. Fzx ^2 ^ I- i^^JT-c? functions 7 /or which there 
exists S{X, ^2) £ TZSX defined in a neighborhood o/tg with the moments 
Hi{t2) for 'y^ computed from i \1.2\} . 

In Section [6] we study the Nevanlinna Pick interpolation problem in the 
class IZSX. Let us recall that in classical Schur analysis, the solution 
of Nevanlinna Pick interpolation problem plays a special role; see for 
instance [FK] , |Dy| , Schur analysis gives a parametrization of all 
solutions (when they exist) for the given data. In the classical case the 
interpolation problem may be formulated in the following way: 

Problem 1.5. Suppose that we are given N pairs of points < wj, sj >, 
j = l,...,Nin]D>x3, where © denotes the open unit disk. Then: 

(1) Give sufficient and necessary conditions, so that there exists 
a function s{z) analytic and contractive in D, and such that 

= Sj, j = 1,...,N. 

(2) Describe the set of all solutions for this problem. 

As is well known, and originates with the work of Pick [P], Problem 
11.51 is solvable if and only if the N x N Hermitian matrix 

1 - SgS* 

Nevanlinna jN] then gave a parametrization of all solutions, i.e., solved 
the second part of Problem 11.51 in the form S{X) = Tq{So{X)), where 
Tq denotes a linear fractional transformation uniquely determined from 
the initial data < Wj, Sj >, and where the parameter S'o(A) runs through 
all functions analytic and contractive in the open unit disk (that is, 
Schur functions). We refer to [BGR] , Dy for more information on this 
classical topic. 

As a generalization of Problem 11.51 in the class IZST we consider the 
following question: 

Problem 1.6 (Nevanlinna-Pick interpolation) . Lei cti, ct2, 7, 7* be ves- 
sel parameters and let I be an interval. Let G N and Wi,i = 1, . . . , N 
be complex numbers. Suppose also that N input functions ^4(^2) satis- 
fying (I2.16P with corresponding spectral parameters Wi 's, and N output 
functions r]i{t2) satisfying (12.171) with the given wt's. 
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(1) Give sufficient and necessary conditions, so that there exists 
S{X,t2) G IZSI such that S{wj,t2)^i{t2) = ^(^2), j = '^,---,N, 
on a sub-interval of I. 

(2) Describe the set of all solutions for this problem. 

One can also consider a harder problem, which can also be considered as 
a generalization of the classical Nevanlinna-Pick interpolation problem, 
which uses the fact that we can also specify the data for different values 
of ts: 

Problem 1.7. Given C^^'^ -valued functions cri,o"2,7 defined on 1, and 
given N quadruples < tljWjy^jjTjj >, where 4 £ -^i ""^i ^ ^j,rij G 
C^^P J = 1,...N, then: 

(1) Give sufficient and necessary conditions, so that there exists 7* 
and S{X, ^2) G IZSX such that S^vujjtDCj = Vj^ j = 1; • • • 1 
on a sub-interval of I containing all the 4 • 

(2) Describe the set of all solutions for this problem. 

If all the values 4 = 4 equal, we have to find a function S'(A, ^2) 
satisfying S{wj, ^2)^3 — Vj- Thus, the above problem is a generalization 
of the classical Nevanlinna-Pick interpolation problem. We also remark 
that we do not address the question of describing the set of all solutions. 

Remarks: The present paper deals with the rational case; the general 
case will be treated in a forthcoming publication. Some of the results 
presented here have been announced in |AM Vj . 

2. ti INVARIANT CONSERVATIVE 2D SYSTEMS. 

The material in this section is taken from [M] and |MVlj . where proofs 
and more details can be found. The origin of this theory can be found 
in the paper [Lij. 

2.1. Definition. An overdetermined conservative ti-invariant 2D sys- 
tem is a linear input-state-output (i/s/o) system, which consists of 
operators depending only on the variable t2 and is of the following 
form: 

' i^x{ti,t2) = Ai{t2)x{ti,t2) + B,{t2)u{t^,t2) 

(2.1) x{tiM) = F{t2,tl)x{ti,tl) + j F{t2,s)B2{s)u{ti,s)ds 

, y{tl-, t2) = U{ti,t2) - B{t2)*x{ti,t2), 

where the variable ti belongs to M, and the variable ^2 belongs to some 
interval I. Furthermore, the input u{ti, ^2) and the output y{ti, ^2) take 
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values in some Hilbert space £ and the state x{ti,t2) takes values in 
Hilbert spaces T-Lt^. We assume that u{ti,t2) and y{ti,t2) are contin- 
uous functions of each variable when the other variable is fixed. The 
operators of the system are supposed to satisfy the following: 

Assumptions 2.1. 



(1) Ai{t2) : 7/(2 — 7- 7it2, and -8(^2) : £ 'Ht2 ^'^s bounded operators 
for all t2, 

(2) The functions di, a2, 7, 7* : £^ — )■ are continuous in the oper- 
ator norm topology. 

(3) o"i(f2) is an invertible operator for every t2 G I. 

(4) F{t, s) is an evolution continuous semi-group. 

For continuous inputs u{ti,t2), the inner state is continuously differ- 
entiable. Requiring now the invariance of the system transition from 
(^5,^2) to (^1,^2) via the points (t?,t2) and (^1,^2), is equivalent to the 
equality of second order partial derivatives of x{ti,t2): 

Substituting in this equality the system equations we obtain that for 
the free evolution u{ti,t2) = the so called Lax equation holds 

(2.3) A,{t2) = F{t2,tl)Mt2)F{tlt2). 

Inserting ( 12. 3p into (12. 2p we see that the input u{ti,t2) has to satisfy 
the following PDE 

B{t2)a2{j:2)i^u{t^,t2) - B{t2)a^{t2)i-^u{ti,t2)- 
{A,{t2)B{t2)a2{t2) + F{t2,tl)i-^[F{tlt2)B{t2)cr,{t2)])u{t,,t2) = 0. 

Assuming the existence of a function 7(^2) satisfying 
(2.4) 

A,{t2)B{t2)a2{t2) + F{t2,tl)-^^[F{tl,t2)B{t2)cri{t2)] = -B{t2h{t2) 
we obtain that it is enough that u{ti, ^2) satisfies the PDE 

(2.5) a2{t2)j-u{ti,t2) - a^it2)—u{t^,t2) + i{t2)u{t^,t2) = 0. 

Oti Ot2 

The output y{ti,t2) should satisfy the output compatibility condition of 
the same type as for the input compatibility condition (12. 5p . namely: 

(2.6) a2{t2)—y{t,,t2) - ^1(^2)7^2/(^1,^2) + 7*(^2)2/(ti, ^2) = 0. 

oti Ot2 
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Inserting here y(ti,t2) = u{ti,t2) — B{t2)*x{t2,t2) we obtain that 

(2.7) = a2{t2)B{t2yA,{t2)F{t2A)- 

-a^{t2)^[B{t2rF{t2A)]+i*{t2mhrF{t2A) 

012 

(2.8) 7(^2) = (r,{t2)B{t2rB{t2)a2it2) - 

-02{t2)B{t2yB{t2)cri{t2)+l,{t2). 

The fact that the system is lossless comes from the requirement of the 
so called energy balance equations: 

d 

— {X{ti,t2),x{ti,t2))nt2 + {Cri{t2)y{ti,t2),y{ti,t2))£ = 

= {<yzit2)u{ti,t2),u{ti,t2))£, 1 = 1,2, 

which means that the energy of the output is distributed between the 
energy of the input and the change of the energy of the state of the 
system. Immediate consequences of this requirement are 

(2.9) = Aiit2) + Al{t2) + B{t2)aiit2)Bit2y, 

(2.10) -^[F*(t2, 4)F{t2, 4)] = F*{t2, tl)B{t2ya2{t2)B{t2)F{t2, t°). 
at2 

In this manner we obtain the notion of conservative vessel in the integral 
form, which is a collection of operators and spaces 

QJ = (Ai(t2), F{t2, t^), B{t2); ai(t2), ^2(^2), 7(^2), 7* (^2); nt„S) 



where the operators satisfy the regularity assumptions 12. and the 
following vessel conditions: 

= A,{t2) + Al{t2) + B{t2ra,{t2)B{t2) O 

||F(t2,t^)x(ti,tO)f -||x(ti,tO)f =^ 

= J^^{a2is)B{s)x{t^,s),B{s)x{t,,s))ds (^M) 

F{t2,tl)A^{t'',) = A^{t2)Fit2,tl) ^i) 

= A.{Fitlt2)B{t2)a,it2)) + 

+F{tlt2)A^{t2)B{t2)a2{t2) + F{tlt2)B{t2)^{t2) (HaD 

= ai(t2)^[5(t2)7(t2,t^)]- 

-(r2{t2)B{t2YAr{t2)F{t2,tl)~^.{t2)B{t2rF{t2,tl) (1221) 

7(^2) = -cy2{t2)B{t2yB{t2)ai{t2) + 

+(Jl{t2)B{t2YB{t2)(J2{t2) + 7*(^2) dH 

In order to simplify some notations we make the following definition 
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Definition 2.2. Let 

U = U(Ai, A2; h) = (T2{t2)Xl - fTi(t2)A2 + 7(^2) 

and similarly 

U* = U*(Ai, A2; t2) = cr2(i2)Ai - o-i(t2)A2 + 7*(^2) 
Then the vessel QJ is naturally associated to the system (12. ip 



S : < 



^^X{h,t2) = Ai{t2) X{h,t2) + B{t2)cri{t2) U{ti,t2) 

x{h,t2) = F{t2,f2Mh,t9,) + J F{t2,s)B{s)a2{s)u{h,s)ds 

4 



^ y{tiM) = u{h,t2) - B{t2y X{ti,t2). 

with inputs and outputs satisfying the compatibility conditions (12. 5 p 
and (12.60 . i.e. satisfy: 

The theory of such vessels, developed in [Ml IMVl] enables to find a 
more convenient form of the vessel. Denoting = T^^o, Ai = ^1(^2), 

F*{t2,4)F{t2,4) =^~\t2) andS(t2) = ^2)^(^2), we shall obtain 
the following notion, first introduced in |M2] . 

Definition 2.3. A (differential) conservative vessel associated to the 
vessel parameters is a collection of operators and spaces 

(2.11) QJ= (Ai,S(t2),X(t2); ai(t2), 0-2(^2), 7(^2), 7*(^2);^,^), 
where the operators satisfy the following vessel conditions: 

(2.12) = ^{B{t2)a,{t2)) + AiB{t2)a2it2) + B{t2h{t2), 

dt2 

(2.13) AiX(t2) + nh)Al = B{t2)a^{t2)B{t2)\ 

(2.14) 4-X(t2) = B{t2)a2{t2)B{t2)\ 

at2 

(2.15) 7*(^2) = 7(^2) + ai{t2)B{t2y^-\t2)B{t2)cy2{t2)- 

-a2it2)B{t2rX-\t2)B{t2)a,{t2) 

This representation of a vessel is the most convenient when one focuses 
on the notion of transfer function, as we do in the next subsection. 
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2.2. Transfer function. Performing a separation of variables as fol- 
lows 

U{ti,t2) =ux{t2)e^'\ 

y{tiM) = yx{t2)e^'\ 

we arrive at the notion of a transfer function. Note that u{ti,t2) and 
y{ti, ^2) satisfy PDEs, but U\{t2) and yx{h) are solutions of LDEs with 
spectral parameter A, 

(2.16) \J{X,^;t2)u{h,t2) = 0, 

012 

(2.17) U(A,^;t2)l/(ti,t2) = 0. 

The corresponding i/s/o system becomes 

f xxit2) = (A/ - A,{t2))-'B{t2)aiit2)uxit2) 

i-xx{t2) = F{t2,4)xxitl) + J F{t2,s)B2{s)ux{s)ds 

^ yx{t2) = Ux{t2) - B{t2yXx{t2) 

The output yxih) = uxih) — B{t2)*xx{t2) may be found from the first 
i/s/o equation: 

y\{t2) = S(\,t2)ux(t2), 
using the transfer function 

S{X,t2) =I-B{t2nXI~A,{t2))-'B{t2)ai{t2) 

= I - S(t2)*(A/ - F{t2, t°) A(tO)F(tO, t2))-'B{t2)a,{t2) 
= I - B{t2yF{t2,tl){XI - A^{tl)y^F{tlt2)B{t2)a^{t2) 

= I - B{t2yF*{tlt2)F*{t2A)F{t2A){^I - A,)-^B{t2)cy,{t2) 
= I - B{t2yX-\t2){XI - A,)-^B{t2)a,{t2) 

and we obtain that 

(2.18) 5(A, t2) = I- B{t2yX-\t2){XI - A^)-'B{t2)a^{t2). 

Proposition 2.4. The transfer function S{X,t2) defined by (12.181) has 
the following properties: 

(1) For allt2, S{X,t2) is an analytic function of X in the neighbor- 
hood of 00, where it satisfies: 

S{oo,t2) = Ip 

(2) For all X, S{X,t2) is a continuous function 0/^2- 
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(3) For A in the domain of analyticity of S{X,t2) : 

(2.19) Six, ta) Vi(t2)^(A, h) < a^ih), > 0, 
and 

(2.20) S{\ t2yai{t2)S{\ t2) = Mt2), 3?A = 0. 

(4) Maps solutions of the input LDE fl2.16p with spectral parameter 
A to the output LDE (12.171) with the same spectral parameter. 

Proof: These properties are easily checked, and follow from the defi- 
nition of S'(A, ^2): 

S{\ t2) = I- B{t2nXI - A,{t2))-'B{t2)a^{t2). 

The function S'(A,t2) is analytic for A > ||yli(t2)|| and since all the 
operators are bounded, we have S{oo,t2) = Ip. The second property 
follows from the regularity assumptions 12.11 The third property follows 
from straightforward calculations: 

5(A,t2)Vi(t2)5(A,t2)-ai(t2) = 

-2^{X)a^{t2)B{t2yCXI - Al{t2))-\XI - A^{t2))-'B{t2)cr^{t2) 

Here the sign of — 3?(A) determines the sign of '5'(A, t2)*c"i(^2)5'(A, ^2) — 
(71(^2) and thus the third property is obtained. The fourth property 
follows directly from our construction. □ 

Remark: When dim'H < 00, we obtain that S{X,t2) is a rational 
function of A for every t2- 

It is an interesting fact that also the converse of Proposition 12.41 holds. 
It is proved in [M], |MVcl chapter 5] . 

Theorem 2.5. For any function of two variables S{X,t2), satisfying 
the conditions of Proposition 2.4\ there is a conservative ti-invariant 



vessel whose transfer function is S{X,t2] 

We define the class of transfer functions mentioned in the introduction 
as follows: 

Definition 2.6 f pVc] ). The class SI = 5X(U(A, — ; t2), U,(A, — ; ^2 

OT2 (yt2 

consists of functions S{X,t2) of two variables, which are 

(1) analytic in a neighborhood of X = 00 for all t2 and where it 
holds S{oo,t2) = Ip, 

(2) continuous as functions of t2 for all X, 

(3) satisfy (12.191) and (12.201) in the domain of analyticity of S , 
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(4) map solutions of the input LDE ^2.10i) with spectral parameter 
A to the output LDE l\2.17\ ) with the same spectral parameter 

Recall (see |CoLej ) that to every LDE can be associated an invertible 
matrix (or operator) function $(^2? ^2)5 called the fundamental solution, 
which takes value / at some preassigned value ^ind such that any 
other solution u{t2) of the LDE, with initial condition u{t2) = uq is of 
the form 

Let <l'(A,t2)^2) ^*(A,t2)^2) be the fundamental solutions of the 
input LDE fl2.16p and the output LDE fl2.17p respectively, where we 
have added in the notation the dependence in A. Then, 

(2.21) Six, t2)$(A, t2, t^) = $.(A, t2, tl)SiX, t°) 

and consequently S{X, ^2) satisfies the following LDE 



(2.22) 



d 



Qt^S{X,t2) = a^\t2){a2{t2)X + 7*(t2))5(A, t2)- 

-S{X,t2)cJ^\t2){cX2{t2)X + ^{t2)). 
3. SCHUR ANALYSIS FOR THE CLASSICAL CASE 



One of the approaches to the tangential Schur algorithm for a rational 
matrix function S{X) G C^^^ is based on the theory of reproducing ker- 
nel Hilbert spaces of the kind introduced by de Branges and Rovnyak; 
see [dBRlj . |dBR2] . Dy for information on these spaces. The paper 



[AD] considers the case of column-valued functions. In this section 
we adapt the results of [AD] to the case of row-valued functions. We 
also present a new formula for a realization of a Schur function after 
implementation of the tangential Schur algorithm. 

3.1. Schur functions and Reproducing Kernel Hilbert spaces. 

In this section ai denotes a fixed self-adjoint and invertible (but not 
necessarily unitary) matrix in C-^^^. Let S{X) be a rational function, 
(Ti-inner in the open right half plane, i.e. 

S{X)*aiS{X) - (Ti < 
at all points in the domain of analyticity i7(S') of S* in C^, and 

S{XyaiS{X) -ai =0 
at all points on the imaginary axis where S is defined. Then, the kernel 

(3.1) Ks{X,w) = —^^ 
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is positive for A, it; G ^2(5'), and the space of rational C^^^- valued func- 
tions 

n 

■t=i 

is finite dimensional. These well-known facts can be proved using re- 
alization theory; see for instance }AGj . Furthermore, 'H(S') is the re- 
producing kernel Hilbert space, associated to the kernel i^5(A, w). The 
inner product is defined by 

{cKs{\,iy),dKs{X,w))ns = {cK s{w , v)) c^xp = cKs{w,u)d*. 
For an arbitrary /(A) G 'H(S') we have the reproducing kernel property 

{f{X),CK{X,w))ns = f{w)C- 

More generally, let now be a finite dimensional Hilbert space of 
C^^'P-valued functions defined in some set fl, and let {/i(A), . . . , /Ar(A)} 
be a basis of Ai. Let X G C^^p denote the Gram matrix with i,j entry 
given by 



(3.2) 



= (fjWi .hW)M, ^,3 = l,•••,p• 



It is easily seen that the space is a reproducing kernel Hilbert space 
with kernel given by the formula 



(3.3) K{X,w)= [ Mw)* 
We set 

(3.4) F(A) 



fNiw)* ] 



/i(A) 



/i(A) 
L fN{X) 



L fp{X) 



Assume now that A4 consists of rational functions, defined on a set 
For a G Q{Ai) the backward-shift operator is defined by 



/(A) - /(«) 

A — a 



Suppose that: 

(1) The space Ai is invariant under the action of R^, 

(2) The functions /i(A) has the property that /i(oo) 
F(oo) = 0, 



0, i.e.. 
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then the function F given by (13 ■4p can be written as 

F(A) = {XI-A)-^Bai, 

for suitably chosen matrices A, B. In this special case, formula (13. 3 p 
takes the form 

(3.5) K{X, w) = (TiB*{wI - A*)-^^-\\I - A)-^Bai. 

As mentioned at the beginning of this section we are interested in 
kernels of the form (13. ip for some ai-inner rational function S. We now 
recall the characterization of these spaces, and first note the following: 
equation (13. ip leads to 

(3.6) ai - S{wYa,S{X) ^ p^^y^-rp^^^_ 

If S is analytic at infinity and satisfies there S'(oo) = Jp, and letting 
w — )■ oo in this equation, we obtain the formula 

(3.7) 5(A) = /p-5*X-i(A/-A)-^Sai. 

Theorem 3.1. Let A4 be a finite dimensional Hilbert space of C^^^- 
valued rational functions, which are zero at infinity. Suppose that 
Ra-M. C Ai for a G and let X be its Gram matrix with re- 

spect to F{\). Then M. = 'H(S') for S defined by (13. 7p if and only if 
the Lyapunov equation 

(3.8) AX + XA* + BaiB* = 
holds. 

When the spectrum of the operator A is in the open left half plane 
3fJA < 0, one has 

where ii2,cri is the Hardy space with the inner product 



We set J 



-ai 
ai 



Note that J is both invertible and self- 



adjoint, and one can define J-inner rational functions. Let G 



On 6x2 

©21 ©22 



be a J-inner rational function; we introduce the linear fractional trans- 
formation 

Te{w) = (On + i^e2i)-'(ei2 + w^e22). 
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Theorem 3.2. Let S and Q he respectively ai-inner J -inner rational 
functions. Then there exists a ai-inner rational function W such that 
S = Tq{W) if and only if the map 



(3.9) 



F 



-S{X) 
L 



is a contraction from 'H(0) to 'H(5'). 

This theorem originates with the work of de Branges and Rovnyak (see 
|dBRlt Theorem 13 p. 305]), where it is proved in a more general 
setting, and is one of the key ingredients to the reproducing kernel 
approach to the Schur algorithm. The proof is the same as for column- 
valued functions, and is omitted. 

As a special case of the previous theorem we have: 

Corollary 3.3. Given S G IZS, Wi,...,Wn G C+, and row vectors 

B= ■■ ■■ , Ai = dmg[-wl,...,-wl]. 

Let M be a Hilbert space of row vectors spanned by the rows of the 
matrix-valued function 

F{X) = {XI-Ai)-^BJ. 
Let X be the solution of the Lyapunov equation 

(3.10) AX + XA* + BJB* = 
and assume X > 0. Let G be defined by 

(3.11) e(A) = hp - B*X-\XIn - Ai)-^BJ. 
Then there exists So E TtS such that S = Te(5'o). 

Proof: Let us denote by t the map ( 13. 9p . For an arbitrary element 
/(A) = r]F{X), we shall obtain that 



tf{X) =r]{XI-A,)-^BJ 



X+Wi 



6 



77 diagi 

<nS{Wn)* in 

E.ATfe-^.(^i-^(«^.)*^i'5(A)) 
Hiiir}iKs{X,Wi) 



(TiS 
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and consequently, 

{tf, tf) = ^ir]iKs{X, Wi),^. CjTiiKsiX, Wj)) = 

= (/,/)• 

Thus t is an isometry and Theorem 13.21 allows to conclude. 



□ 



The following result shows that the assumption X > in the statement 
of Corollary 13.31 can always be achieved for ri = 1. 

Lemma 3.4. Given S G TZS which is not the function identically 
equal to Ip. Then there exist a pair {^,w) G C^^^ x C+ such that the 
corresponding X > 0. 

Proof: We proceed by contradiction. Assume that for each w G C+ 
and for each vector ^, it holds that 

or, equivalently, 

CKs{w,w)C = 0. 

Then, ^^^(A,^;) = 0, and for each / G HiS) 

^f{w)=<f{X,^Ks{X,w) >=0. 

The space 'H(>S') is thus trivial, and its kernel is zero: 

a, - S{wya,S{X) _ 
X + w 

from where we conclude that for each w, X 

S{X) = a^^S-*{w)a^. 
and consequently, S{X) = Ip. □ 

3.2. Analysis of the tangential Schur algorithm. For n = 1, the 

matrix function in (13. lip becomes 
(3.12) 
0(A) 



hp- B*X-\XI - A)-^BJ 



h, 



2p 



-ax 

ai J X(A + w* 



X(A + w* 
"x(A + ty*) 



X(A + w*) 



X(A + u;*) . 



en 912 

©21 ©22 



ON THE CLASS nSI 

where using the Lyapunov equation (13.101) 

+ Xv4* = X(-u;* -w) = 

= -BJB* = - [ -?7 e ] 
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" -(Tl " 




—7]* 


ai 







we find that 
(3.13) 



X 



^(jiC - V(JiV* 
w + w* 



To the best of our knowledge, Theorem 13.51 below is new. It will be of 
much use in the following sections. 

Theorem 3.5. Let Sq G IZS with minimal realization 

(3.14) So{X) = Ip- B*X,\XI - Ao)-'Boa^, 

and let Q be given by ( 13. 12^ . Then S = Tq{So) G IZS and a minimal 
realization of S is given by 



SiX) = Ip- B*sX-s\XI - As)-'Bsai, 



where 
(3.15) 

(3.16) 



B, 



Xc 



Bo 

Xo 
X 



(3.17) 



A, 



Ao 

-r]aiB*Xo' -w* 
Proof: From the definition 



Xo 

w + w* 
X, 

X 



S{X) =Te{So{X)) = 

= (Gn + So{X)Q2i)-\Qi2 + So{X)Q22) 



r] rjcri 



SoiX): 



X 



X{X + uj*) 'X(A + «;*)' 

X(A + ^*) ^"^"^^^^^^ X(A + ^*) 
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Let US denote here So{\)^* = a*; then the preceding expression be- 
comes 

r]*riai a*rj(Ti 



S{\) = (/,+ 



X(A + w*) X{\ + w*)' 



[t]* - a*)ria 



1^-1 



X(A + w*)' 



Using (/ + 6*a)-i = I-b*a 
S{X) = (/,- 



X(A + i(;*) 
1 



(^o(A) 



X(A + w* 



1 + a6* 
(?7* - a*)?7cri 1 



we get 



X(A + «;* 



770-1(77* - a* 



-)(^o(A) + 



(77* - a*)^eri. 
X(A + «;*) ■ 



5'o(A) 



X(A + m;*) 



X(A + w*) + 77^1(77* - 5'o(A)e) ' 

Suppose further that there is a reahzation of the form (I3.14p for the 
function 5*0 (A). Inserting it here we shall obtain 

c(\\ ^ c (\\ ~ 'S'o(A)r)[Wi^o(A) - g^i] _ 

^ ' °^ ^ X(A + 7/;*) + 77^1(77* -S'o(A)e) 

= Ip-B*Xo\XI-Ao)-'Boai- 

(77*-^o(A)r)[wi^o(A)-e^i] 

X(A + w*) + 77^1(77* - So{X)^*) ■ 

Let us denote 

M =X{X + w*) + r]ai{7]*~So{X)e) = 

= X(A + w*) + 77^1(77* - C) + B*oX,\XI - Ao)-'Boa,C- 

The preceding formula for S'(A) becomes 

(77*-5o(A)r)[wi^o(A)-e^i] 



SiX) = Ip-B*X,\XI-Ao)-'Boa, 



Xn^ 



^0 





X(A + w* 

X 



1 



-5a~^ 1 



77^1(77* -5o(A)e) 



Bo 
77-^ 



where 

a = XI — Aq 

[5 =-Boa,e 
6 =r]a^B*X^\ 

By definition of M we have 

M = X(A + w*) + 77(71(77* - r) - 5a-^/3, 
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and hence we obtain from the formula 
H -1 



a P 
5 D 



1 



D - 6a-^/3 
that the last expression for S{X) is 



-6a~^ 1 



S{X) 



1 



X 



X,' 

M L ■' ^ J [ 1 

-5a~^ 1 



X 



X, 



-1 











X 



a 



/5 





-1 











A/ -A 



X 






X 



0"! 



X 



X 



-1 



B, 



ai. 



Thus we have obtained the realization fl3.15p - fl3.17p . Furthermore, the 
Lyapunov equation fl3.10p holds since 



AsXs + XsA*s 

An U 

+ 



X 



Xo 
X 



+ 



Xo 
X 



X 



X ^ _ 



X 

AqXo + XoA* -Bo(Ti{r]* -C) 



-Bsa.B^ 



s- 



The last equality follows easily from the Lyapunov equation for the 
given realization of and the formula fl3.13p . □ 
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4. The tangential Schur algorithm in the class IZSI 

Our strategy to the tangential Schur algorithm in the class IZSI re- 
lies on the following theorem. This theorem shows in particular that 
one cannot use the naive approach of applying the classical tangential 
Schur algorithm for each t2 (that is, looking at t2 as a mere parameter). 
The theorem itself is stated and proved for the non-rational case. 

Theorem 4.1. Let us fix the parameters <Ji,a2, and 7, and the interval 
I. Then for every t^ & I there is a one-to-one correspondence between 
pairs (7*, S) such that S G SI and 7* continuous in a neighborhood of 
^2, and functions Y{X) G «S(t2)- 

Proof: Let (p be the map which to a pair (7=,,, S) associates the func- 
tion S{X,t2) G «S(t2)- The converse map ip '■ ~^ y{^,h) was 
introduced in [M], |MVll chapter 7] in a more general setting, and is 
defined as follows. Suppose that we have realized the transfer function 
Y[X) in the form 

y(A) = Ip- B*X,\XI - A,)-'Boa,{tl). 

Then we construct -8(^2) from the differential equation with the spec- 
tral matrix parameter Ai: 

(4.1) ■^[B{t2)a^{t2)]+A,B{t2)a2{t2)+B{t2Mt2) = 0, B{tl) = B^. 
0x2 

In fact, the function ^(^2) is given by the formula 

(4.2) B{t2) = j{XI - A,)-'Boa^^^-\X,t2A)d\. 

Next we construct X()f:2) on the maximal interval X, where it is invert- 
ible (12. lip via the formula 

(4.3) 4-X(t2) = B{t2)a2{t2)B{t2)\ X{tl) = Xq. 
dt2 

Finally, we define 

7*(^2) = 7(^2) + (y2{t2)B{t2YX-\t2)B{t2)ai{t2)- 

-a,{t2)B{t2)*X-\t2)B{t2)a2{t2). 
Then easy computations show that the function 

(4.4) S{X,t2) = Ip- B{t2rX-\t2){XI - Ai)-'B{t2)cr,{t2) 

is in the class CX corresponding to the parameters cri(t2),0"2(^2),7(^2), 
and 7=1,(^2) for t2 G I. Moreover, the Lyapunov equation (13. 8p holds for 
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every ^2 £ I, and thus the reahzation (14.41) is minimal for every it2 £ I- 

Note that the composition (p o tfj = id, since starting from a function 
Y G S{t2), constructing Y{X,t2) and taking its value at we shall 
obtain again Y from the initial conditions of the differential equations 
(CT and (gJl) defining 5(t2),X(t2). 

In order to show that ifjocj) = id, we start from a function S{X, ^2) £ «5X 
and take its value S{X,t2) G ^{^2)- Using the construction above, we 
shall obtain a function Y{X,t2)- Note that the two functions S{X,t2) 
and y (A, ^2) have the same value at ^2 and maps solutions of the same 
input LDE to (possibly different) output LDEs, i.e. : 

SiX,t2) = $,(A,t2,t^)^(A,t°)$-HA,t2,t^), 
Y{X, t2) = $:(A, t2, 4)S{X, t°)$-i(A, t2, tl). 

Then the function S~^{X,t2)Y{X,t2) is equal to Ip at infinity and is 
entire. By Liouville's theorem it is a constant function and is equal to 
Ip. Thus 

^,{x,t2,4) = ^:{x,t2,4), 

from where we obtain that 

^:\x,t2,4)Ki>^,t2,ti) = ip. 

Differentiating both sides of this last equation we get to 

=^^[K\x,t2,4mx,t2,m = 

= <l>~\X,t2,t''2)c^^\t2)i-7{t2)+iit2)mX,t2,4)- 

Since the matrices $=„(A, ^2, ^2)1 ^2, ^2)5 '^1(^2) are invertible we ob- 
tain that 7(^2) = I'ih)- □ 

Remark: Last theorems claims that the correspondence is between 
"initial" values ^(Ajtg) and pairs (S'(A, ^2), 7*(^2))- Notice that it is 
possible to obtain functions with the same 7*(t2) with different initial 
conditions: 

Proposition 4.2. Suppose that there exists a function Y G S{t2), 
which commutes with $(A,t2,i2) ^'^^ suppose that a function S G «SX 
corresponds to certain vessel parameters (T2, 7, 7*. Then the func- 
tion SY belongs to the class SX and corresponds to the same vessel 
parameters ai, 02, 7, 7*- 

Proof: Using formula 12.211 we obtain that 

S{X, t2) = $.(A, t2, t^SiX, tl)^-\X, t2, t'2). 
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Consequently, 

SiX,t2)YiX) = ^4X,t2,4)SiX,4)YiX)^-\X,h,i) 

intertwines solutions of the input fl2.16p and the output (12.171) ODEs 
with the spectral parameter A, and is identity at infinity, because 5" 
and Y and their product are such. Thus by the definition the function 
SY G ST and corresponds to the same spectral parameters as S. □ 
For example (to be studied in subsection 15. 2p . taking 





"01" 




"10" 




"00" 


dl = 


1 


,0-2 = 





,7 = 


i 



we shall obtain that 



^x,h,tl) = v- 



g-A;(t2 





-4) 





=fc{t2-t§) 



where 



V 



V2k 
X~ 



i)Vx 



Taking Y, which commutes with ^{cr2X + 7) 
form 




A 



i.e. of the 



Y{X) 



a(A) 
c(A) 



ic{X) 
a{X) 



we shall obtain that for any S G SI, the function SY G «SX and 
corresponds to the same vessel parameters. 

A stronger condition on 7*(t2) will be presented in sectional after we 
shall study Markov moments of the transfer functions. 
We now focus on the rational case. We first note that starting from dif- 
ferent realizations at we shall obtain the same 7*(t2)- More precisely 
we have the following theorem: 

Theorem 4.3. Suppose that there are two minimal realizations of the 
function S{X,t2) G IZS: 

S{X,tl) = 1,- B*XJ\XI - A^)-'B,a,{tl), £ = 1,2, 

with associated similarity matrix V . Then the functions 



SeiX, t2) = Ip- Beit2rXi\t2)iXI - Ae)-'Beit2)a,it2), 



1,2, 



obtained from these realizations via the construction in Theorem \4-l 
are the same, and this holds if and only if: 

(4.5) B2{t2) = VB,{t2), 
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where t2 varies in a neighborhood of t^- Moreover, in this case the 
following formula holds 

X2(t2) = VX^{t2)V*. 

Proof: Equality of the two reahzations means that there exists an 
invertible matrix V such that 



^2 1 



A2 = VAiV-\ B2ai{tl) = B^X^^V = B^X^ 

from which follows (see |AGt Lemma 2.1 p. 184] for instance) that 

X2 = VXiV*. 

By the construction described in Theorem 14. 1^ the function Bi{t2) will 
satisfy (gl]) 

-^[Si(t2)ai(t2)] + AiBi{t2)a2it2) + 51(^2)7(^2) = 0, Bi{tl) = Bi, 
0X2 

and the function -82(^2) will satisfy the same equation with A2 instead 
of Ai: 

-^[B2{t2)a^{t2)] + A2B2{t2)a2{t2) + B2{t2)l{t2) = 0, B2{tl) = B2. 

at2 

Using the equahties A2 = VAiV^^, B2 = VBi we obtain that the 
function V~^B2(t2) satisfies 

^^[V-'B2{t2)ai{t2)] + AiV-'B2{t2)a2{t2) + ^-152(^2)7(^2) = 0, 
V-'B2it^2) = Bu 

which is the same differential equation as for Bi{t2). Thus ^2(^2) = 
VBi(t2). Similarly, considering the differential equations 

-^X,(t2) = B,{t2)aMt2)\ X,{tl) = X„ z = 1, 2, 
we obtain that 

X2(t2) = VX^{t2)V*. 

Consequently, 

B*2{t2)X2\t2)B{t2) = Bl{t2)V*V-'*X^\t2)V-'VB,{t2) 
= Bl{t2)X^\t2)B^{t2), 

from where we conclude that the same function 7* is associated to 
Si{X,t2) and S'2(A,t2)- Since these functions coincide for t2 = ^2 
map the same input ODE, we obtain that they are equal in a neigh- 
borhood of t2. □ 

The following notion has been introduced and studied in |M2j . 
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Definition 4.4. Let S G IZSX with a realization MAh . The function 

T{t2) = detX(t2). 
is called the r-function associated to S . 

It follows from Theorem 14.31 that the r function is well defined up to a 
multiphcative strictly positive constant. Indeed using the notation of 
the theorem, 

detX2(t2) = det[VXi{t2)V*] = detXi{t2) det{VV*). 

We now introduce the counterpart of the tangential Schur algorithm 
in the class IZSI. Let S G IZSI and fix ^ I- The function 
^(A, ^ 7?.«S. Consider now a space M with X > and corresponding 
function O as in Corollary l3.3[ This is always possible in view of Lemma 
13. 4[ It follows from Theorem 13.21 that there exists 5*0 G IZS such that 

(4.6) ^(A,t°)=Te(A)(5o(A)). 

Applying Theorem 14.11 to 5*0, we obtain a uniquely defined function 
'S'o(A,t2) £ TZSI, such that at the relation (14. 6 p holds. 

Definition 4.5. The map T^ q 

S{X,t2)^So{X,t2) 

is the time-varying counterpart of the linear fractional transformation 
(14.61) . We will call it a generalized linear fractional transformation. 

If S'(A,t2) € <5X corresponds to the vessel parameters ai, (T2, 7, 7*, 
then So{X,t2) G «SX corresponds to the vessel parameters cti, a2, 7, 
7o,* for a uniquely defined function 70,* (^2)- Moreover, for ^2 = ^2 
have the usual linear fractional transformation (14.60 . 

As a consequence of Theorem 14.11 the following lemma holds. 

Lemma 4.6. For a given J -inner function and a given point t2, the 
map T^ fO is one-to-one from IZS into IZSI. 

Proof: Notice that for a given t°2 the map Te is injective. Further- 
more, using Theorem 14.11 every S G IZSI is uniquely defined by the 
function 5'(A, fg) ^ IZS. The result follows. □ 

Suppose now that we start from 5*0 (A) = Ip and apply n linear frac- 
tional transformations for a fixed t2, using the data < Wi,^i,r]i > 
{i = l,...n) to construct the corresponding J- inner functions. We 
obtain a function 

Sn{X) = Te(/) = / - B:X-\XI - A^y'B^a, G TLS. 
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Using iteratively formulas fl3.15p . f l3.16p . f l3.17p we obtain that 

Bn = : 
X„ = diag[Xi, . . . ,X„], 
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and 

(4.7) 



m^ii-nl -il) 



r?20-i(r)* - c*) 



iV2 - 



-W„ = 



where Xj is defined by (13.13P : 



X,- 



Wi + w* 



1, . . . Ti. 



Assume now that, starting from the identity matrix we apply this pro- 
cedure for two different sets of data (with the same n) at two different 
points tl and t^- The following theorem answers the question as when 
we obtain the same function, that is, when do we have: 



Theorem 4.7. Suppose that there are given two sets of n triples 
< w},^l,r]l > and < Wi,^f,rif >, with corresponding 6^,£ = 1,2. 
Then necessary and sufficient conditions for equality of the two func- 
tions 



1,2 



are: 



(1) The corresponding matrices A\ and A\ defined by (14.71) are sim- 
ilar, i.e. there exists an invertible matrix V such that A]^ = 
VAlV-\ 

(2) §{\I - A\)-'Bia-,'^-\\tlA)dX = VBl. 



Proof: From Theorem 14. H a necessary and sufficient condition for the 
functions to be equal is that 

s,{\tl) = S2{\A)- 

From 14.31 this holds if and only if 

Al = VAiV-\ Blitl) = VBl Xlitl) = VXlV* 
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for a uniquely defined invertible matrix V. The result follows using 
formula (jtg). □ 

A more general construction in this setting is obtained if one supposes 
that at each step different values of ^2 are chosen. In this case the 
construction of the function S'„(A,t2) is more complicated, and can be 
computed recursively. The formulas are very involved in this case and 
we can see no real advantage to develop them at this point. 

4.1. Linear Fractional Transformations in terms of intertwin- 
ing positive pairs. Suppose that we are given a data of the NP in- 
terpolation problem 11.61 Following the notations of corollary 13.31 let us 
write 



Let us denote by capital Greek letters the following vessel parameters 



Sl(t2 



S2(t2 



-(ri{t2) 
ai(t2) 

fT2(t2) 

^2(^2) 

' l*{t2) 
7(^2) 



J, 



then simple calculations show that B{t2) satisfies (14. ip 

-^[S(t2)Sl(t2)] + A5(t2)S2(t2) + B{t2)T{t2) = 0, = B. 

dt2 

Suppose that X(t2) > is a solution of 

A,X{t2) + Mt2)Al + B{t2)^i{t2)B*{t2) = 0, 



_d_ 

dti 



X{t2) = B{t2)^2{t2)B*{t2), 



which is always possible if ^ for each i = 1, 
following collection 



, n. Then the 



Qj = {A,s(t2),x(t2);Si(t2),S2(t2),r(t2),r,(t2);C^";^^©^}. 

is a vessel for r^,(t2) defined from the linkage condition fl2.15p 

r,(t2) = r(t2) + Si(t2)S(t2)*X-l(t2)S(t2)S2(t2) " 
-S2(t2)5(t2)*X-l(t2)5(t2)Si(t2). 
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Transfer function of the vessel QJ is 

e(A, h) = hp - B{t2rX-\t2){XIn - Ai)-lS(t2)Si(t2), 

which is in 72.<SX(U* © U, U) for 

U = AS2(t2)-Si(t2):^ + r,(t2). 

dt2 

If we denote further the decomposition of 6(^2) as 

e(A,t2) = 
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eii(A,t2) ei2(A,t2) 
e2i(A,t2) e22(A,t2) 



then if one defines 5*0 (A, ^2) so that 

(0n + 5'oe2i)-^(ei2 + 5-0922) = 5, 

then the function ^(A, usually does not intertwine solutions of LDEs 
with spectral parameter A. 
Instead, we define 



W{X,t2)= [ W,iX,t2) W2iX,t2) ] 



SO that 



WiQn + W2Q21 



W1Q12 + ^^2622 = S, 



then the following lemma holds 



Lemma 4.8. The pair of functions W{X^t2) is in «SX(U,U*) and 

— ^ ' 2) i( 2) — {^^^ 2) > g domain of analyticity of W{X,t2) ■ 

A + /i 

Proof: Since 0(A,t2) is invertible for all A out of the spectrum of Ai, 

yx{t2) 



U\{t2 



an element of U is of the form 0(A, 12) 
satisfy fl2.17p and (12.161) respectively. Then 

W{XM)Q{XM 

= [ Wi{XM) W2{XM) ] 

= [ WiOn + W2Q21 mei2 + 1^2622 ] 



where yxit2), ux{t2) 



eii(A,t2) ei2(A,t2) 

e2i(A,t2) e22(A,t2) 



yx{t2) 

U\{t2) 

yx{t2) 
yx{t2) 

U\{t2) 



[ I S{XM) ] 



yx{t2) 

U\{t2) 



yx{t2) + S{X,t2)ux{t2) e U,, 



since yx{t2) and S{X,t2)ux{t2) are in U^,. 
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From the formula WQ = [ ^ S ~\, it follows that 

w{\,h)= [ I s(\.h) 

Consequently, the e.p.e=.on m\t,r^.(t.)m,M ^^^^.^^^^^ 

A + /i 

the domain of analyticity of W{-,t2) becomes 

W{x,t2)j:i{t2)w*{fi,t2) 



[ I SiX,t2) ] 



e-\x,t2)J:l{t2)Q-'*{^^,t2) 

X + fM 



I 

S*{fl,t2) 



Since G(A,t2) is a transfer function of a conservative vessel, its inverse 
is a transfer function too and satisfies 



e-i(A,t2)Si(t2)e(//,t2: 



-1* 



A + ^ 



> Si(t2) 



and consequently, 

W{X,t2)^^{t2)W*{fi,t2) 
X + ft 



> [ I S{X,t2) ]Si(t2) 

> [ / SiX,t2) ] 



I 

S*{fi,t2) 

-a,{t2) 
ai(t2) 
> S{X,t2)ai{t2)S*{fi,t2) - ai{t2) 



I 



1*1 



-x,u 



S*i-X,t2) 



Notice that 

W{x,t2)j:i{t2)w{-x,t2) = 
= [ I s{x,t2) ]e-HA,t2)Si(t2)e- 

= [/ 5(A,t2) ]Si(t2) S*{-X,t2) 

= S{X,t2)ai{t2)S*{-X,t2) - a.ih) = 0, 

by the properties of transfer functions for vessels. 
As a consequence of this theorem, we define 

Definition 4.9. A pair of functions 

W{X,t2)= [ Wi{X,t2) W2{X,t2) ] 

is called positive if the conditions of lemma \4-8\ hold: 

W{X, t2) e 5X(U, U,), W{X, t2)JW{X, t2) > on C_ 



□ 
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5. Markov moments and partial realization problem in 

the class 'ksi 

Let S G TtS, and consider the Laurent expansion at infinity 

oo ^ 

S{\) = Ip- B*X-\XI - A)-^Bai = - ^ _.B*X-^A'Bai 

i=0 

in terms of a given minimal realization. The matrices Hi = B^H'^A^Bcfi 
are called the Markov moments of S. The Partial Realization Problem 
(or moment problem at infinity) in TZS is defined as follows: Given 
the first n + 1 Markov moments Hq, . . . , Hn, find all functions (if any) 
S G IZS with these first n + 1 moments. See |GKLj for a general study 
of the partial realization problem. Similarly, one can define the Markov 
moments for an element S G IZSX. We now give necessary conditions 
which the moments of a function S G IZSI have to satisfy. Then we 
consider the partial realization Problem II. 4[ that is, given n + 1 C^^P- 
valued functions Ho(t2), . . . , i^„(t2), for a fixed choice of cti, a2, find all 
S G IZSI with these first n + 1 moments. 

It is important to notice the following: fixing t2 = solving the 

corresponding classical moment problem will not lead to a solution of 
the problem in the class IZSI because we cannot obtain the function 
7 (and hence 7*) from this solution. 



5.1. Restrictions on Markov moments for functions in TZSI. 

We study the Markov moments of a function S G IZSI, which maps 
solutions of the input ODE f l2.16p to the output ODE fl2.17p using a 
minimal realization (14. 4p of S. First we develop an analogue of the 
formula for 2D vessels with constant coefficients. 

Theorem 5.1. For fixed ai,a2,j, a necessary condition on 7* so that 
there exists a vessel with the vessel parameters cTi,cr2,7,7* is 

det $,(A, t2, tl) = det $(A, ^2,^2) 

Proof: Let S G SI be a function corresponding to the parameters 
'71,0-2,7,7*. Then using a realization (14. 4p 



S{X, t2) = Ip- B{t2rX-\t2){XI - A,)-'B{t2)a,{t2) 
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and Lyapunov equation (13. 8p we shall obtain that 

det S{X,t2) = det (ip - S(t2)*X-i(t2)(A/ - Ai)-iS(t2)ai(t2)) 
= det (/ - S(t2)ai(t2)5(t2)*X-i(t2)(A/ - 
= det J + (AX(t2) + X(t2)A*)X-i(t2)(AJ - 
= det (J + Ai{XI - + X(t2)^tX-i(t2)(A/ - A^)-^) 
= det ((A/ + X{t2)AlX'\t2)){XI - v4i)-i) 
= det(AJ + Al) det(A/ - A)"^ 
= det5(A,t°). 

Consequently, taking determinant of the formula fl2.21l) 

si\,t2Mx.,t2,4) = ^.ix,h,4)six,4) 

we shall obtain that det $^,(A, ^2, ^2) — det $(A, ^2) ^2) points A, 

where det S{X, ^2) exists and is different from zero. Since it happens for 
all points outside the spectrum of Ai and the functions det $*(A, t2, ^2)) 
det $(A, ^2, ^2) entire they are equal for all A. □ 
Since Ai is a constant matrix, at some stage the elements I , Ai, . . . , A^ 
will be linearly dependent and we obtain: 

Lemma 5.2. Given S G IZSI, with Markov moments Hi{t2), i = 
0,1,.... Then exists N and constants jij , j = 1, . . . , N such that 

N+l 

(5.1) Yl ^'jHn^M = 0, n > iV + 1. 

j=0 

We notice that the first moment i^o(^2) satisfies the linkage condition 
(O): 



7*(i2) - 7(^2) = 0-2(^2)^^0(^2) - cri{t2)Ho{t2)a^\t2)a2{t2). 
Let us denote 

^^0(^2) = Cit2)B{t2)a,{t2) = {B{t2)YX-\t2)B{t2)ar{t2). 
The functions B{t2),C{t2) satisfy the following differential equations 

(5.2) -^[B{t2)a^{t2)] + A^B{t2)a2it2) + ^(^2)7(^2) = 

0,12 

(5.3) a,it2)^Cit2) - a2(t2)C(t2)A - 7*(^2)C(f2) = 0, 

0(12 

see [Ml IM Vl] . Thus, differentiating i^fo(^2), we obtain 

^i/o(t2) = ^^[C{t2)Bit2)a,{t2)] 

= cr^\t2)a2{t2)C{t2)AiB{t2)a,it2) - C(t2)AlS(^2)a2(t2)^ 
+ai\t2h,{t2)Ho{t2) - Ho{t2)a^\t2hit2). 
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In other words the second moment Hi(t2) = C{t2)AiB{t2)ai{t2) satis- 
fies the following differential equation 

a^\t2)a2{t2)H,{t2) - Hi{t2)a^\t2)a2{t2) = 

= i-^Ho{t2) - arS*^^o(t2) + Ho{t2)a^\t2)i{t2). 

In the same manner the moments Hi(t2) and i/j_|_i(t2) are connected 
by the following differential equation 

a^\t2)a2{t2)H,+i - H,+i{t2)a^\t2)a2{t2) = 

^^•^^ = ^i/.(t2) - (yi\t2)-i.{t2)H,{t2) + H.,{t2W{t2)l{t2). 

Notice also that 

Hia^^H*ai = B*X-^A{BaiB*X-^Bai = 

= B*X-^A\[-AiX - XAl]X-^Bai = 

= -Hi+i - B*X-^A{XAlX-^Bai = 

= -Hi+i - B*X-^A\-^[-BaiB* - XAl]AlX'^Bai = 

= -Hi+i + B*X~^^^^BaiB*A\X-^Bai + 

+B*X-^^^^XA\A\X-^Bai = 

= -Hi+i + Hi^ia^^Hlai + B^X'^^^^XAfX'^Bai = 

= -H,+i + H,_ia^'H*ai - Hi_2a{^H2ai+ 
+ --- + (-l)Vr'^^m^i 
Consequently, the following formula holds: 

i 

(5.5) i/.+iar^ + i-iya^'H*^, = J^i-^Y^'H^^-.^i^H; . 

j=0 

Finally, we show how the third condition in Proposition 12.41 is reflected 
in the moments Hi{t2)- This condition means that the function S{X, ^2) 
is a"i(t2) contractive, and, for example, the flrst moment -^^0(^2) satisfles 
Ho(t2)cri^(t2) > 0. Using the minimahty property, we have 

spanA'^BCP = C^, 

where P is the dimension of the state space. Thus, in order to have 
X > (or X^^ > 0) it is enough to demand 

5*(A*)"X-U"5 > 0, for each n. 

Using the Lyapunov (13. 8p equation iteratively we shall obtain that this 
condition becomes 

Hoa^^ > n = 0, 

-H2a^^ - a^^H*aiHia^^ > n = 1, 

H^a^^ + a^^H*aiH^a^^ - a^^HlaiH2(J^^ > n = 2, 
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More generally, for each n 

n-l 

(5.6) (-i)"//2n^ri + ^(-i)^+vrii/;aij/2n^i-^ar' > o. 

Theorem 5.3. Suppose that we are given moments Hi{t2), defined in 
a neighborhood of the point ^2 ^ I. Then there exists Hq < such that 
each element Hi^i is uniquely determined from HQ,...,Hi using the 
algebraic formulas flS.Sp . andn^ LDEs with arbitrary initial conditions, 
obtained from (15.41) . 

Remarks: 1. This theorem is of local nature. The number "no" ap- 
pearing in the theorem may vary with the point ^2; but is unchanged 
in a small neighborhood of by continuity. 

2. We emphasize that in order to generate moments one needs produce 
no initial conditions for each moment Hi. In other words, one need 2no 
complex numbers to determine Hq,Hi, etc. 

3. The proof of the theorem allows to produce an algorithm to deter- 
mine explicitly, up to Uq initial conditions, the Markov moments. The 
arguments in the proof of the theorem are illustrated on an example 
in the following subsection. This example exhibits all the difficulties 
present in the general case. 

Proof of I5.3t From formula (15. 5p we obtain that the real or imaginary 
part of i^i+i(t2) is uniquely determined from all the previous moments. 
Suppose that 3fJi/j+i(t2) is known (that is, i is even); then using formula 
(15.41) in which we insert the formula for ^Hi^i{t2) we shall obtain 
an algebraic equation which gives p^ — no relations on the elements 
of QHi^i{t2), where no is the dimension of the kernel of the linear 
operator defined by the left handside of (15.41) . Let us demonstrate it 
more explicitly, let us rewrite the equation (15.41) as a system of p"^ linear 
equations in variables '^H^l^, k,l = 1, . . . ,p'^: 

(5.7) < . 

Here a^i, jS^^i, are functions of t2 derived from the vessel parameters. 
Using basic algebra manipulations, it is enough to find a maximally 
independent subset of equations in the left side of (15. 7p . and to obtain 
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a system of — uq independent equations 
(5.8) 

for linear functions fj and a system of linear dependent ones 

at2 

0,12 

= 9no{^'^m2),'^H,{t2),m,+^) 
dt2 

for linear functions Qj. 

Consequently, from flS.Sp we express — no elements of Q'i/j+i as func- 
tions of other tiq elements of -ffj+i and of Hi. On the other hand, using 
the same formula (15.41) with z + 1 plugged instead of z, we shall obtain no 
differential equations for elements of which can be considered 

as restrictions on no unknown elements: 

(5.9) = ^7,(-^$>ff,+i(t2), 53//,+i(t2)), J = 1, . . . , no. 

Notice that there are at least no unknown elements between these equa- 
tions, because of the appearance of -— 53ifi+i(f2) at the functions Qj. 

dt2 

Since we obtain algebraic (15. 8 p and differential (15. 9p equations, they 
are all independent. So, if all elements of QHi^i appear at these equa- 
tions we shall obtain that one can uniquely solve them up to no initial 
conditions. 

If this is not the case, some of the elements, say Pq < hq will not appear 
in the equations (15. 8p and (15. 9p . We sum up these consideration with 
the following description: we have obtained that p"^ — po elements of 
'^Hi^i are determined from po unknown elements of Si/j+i and from 
Hi via p"^ equations (15. 8p and (15. 9p . 

By induction the same will hold for -ffi+2- In other words, — po 
elements of Hi^2 are determined from po unknown elements of i/i+2 
and from i/j+i via p"^ equations (15. 8p and (15. 9p . This means that po 
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relations are obtained for the elements of which are actually 

equations on the elements of -ffj+i- This produce differential equations 
on pq unknown elements of i^i+i- Notice that the unknown elements 
will appear only with one differentiation, since their derivatives did not 
appear at the previous stage. Thus i^j+i is uniquely determined up to 
no initial conditions for some of its elements. □ 
Next theorem put some light on the connection between equations, 
which determine the transfer function S{X, ^2) and its moments. This 
theorem is similar to the property of a solution of a Riccati equation 
[Zel theorem 2.1] 

Theorem 5.4. Suppose that S{X,t2) is a continuous function 0/^2 
for each X, meromorphic in X for each ^2 and satisfies S{oo,t2) = I- 
Suppose also that S{X, t2) is an intertwining function of LDEs fl2.16l) 
and fl2.17p . Then if the equality 

S{XM) = (yi\h)S-^*{-XM)^iit2) 
holds for t\, then it holds for all t2- 

Proof: Since S{X,t2) intertwines solutions of (12.161) and (12.171) . then 
it satisfies the differential equation (12.221) 

-^SiX,t2) = a^\t2){a2it2)X + ^,{t2))SiX,t2) - 
012 

-S{X, t2)a^\t2){a2{t2)X + 7(t2)). 

Consequently, using properties of 7*, 7 appearing in definition 11.11 we 
obtain that the function erf ^(t2)5'~^*(— A, t2)o'i(^2) satisfies the same 
differential equation. If these two functions are equal at from the 
uniqueness of solution for a differential equation with continuous coef- 
ficients, they are also equal for all t2- □ 

Corollary 5.5. Suppose that S{X,t2) satisfies conditions of theorem 
5.4\ then the moments the equation for the moments (15.51) is excessive. 

Proof: From theorem 15.41 it follows that 

SiX,t2)a^\t2)S*i-X,t2)a,it2) = I, 

where taking the expansion in moments for S{X,t2) and for S*{—X,t2) 
we will get the formulas (15. 5p . □ 
We want to present next a necessary restriction on 7*(t2)j derived from 
the existence of a finite dimensional vessel: 

Theorem 5.6. Let cri,cr2,7,7* be vessel parameters, and QJ a finite 
dimensional vessel corresponding to the m parameters. Then the en- 
treis of the function 7* satisfies a polynomial differential equation of 
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finite order with coefficients in the differential ring TZ, generated by 
(the entries of) ai,a^'^,a2,'j. 

Proof: Suppose that the transfer function of the vessel 03, defined in 
dMH) is 

oo 

S{\, t2) = I- B*{t2)X-\t2){XI - A,)-'B{t2)a,{t,) = J - ^ 

i=0 

and the hnkage condition is 

(Ti\hh,{t2) = cj^\t2h{t2) + [cr^\t2)a2{t2),Ho{t2)]. 

Notice that if we differentiate this formula and use the equation for the 
derivative of i/o(^2) from the equation (15.41) . we shall get 

^K'(^2)7*(i2)) = ^K'(t2)7(t2)) + [^^icT^\t2)a2it2)),Hoit2)] + 

+ [a^\t2)a2{t2),^H^{t2)] = 

= ^{cr^\t2Ht2)) + foo{Ho{t2)) + foi{Hi{t2)) 
ai2 

for linear in moments functions /oo, foi with coefficients depending on 
IZ and 7=K. Similarly, differentiating this expression and using formula 

(15.41) for - — i^o(^2) and - — Hi{t2), we shall obtain that the second 

Ut2 dt2 

derivative is 

^(ar'(t2)7*(i2)) = 

= ^K'(t2)7(t2)) + flo{Ho) + fn{H,{t2)) + fl2{H2{t2)). 

for linear in the moments functions /io,/ii,/i2 with coefficients de- 
pending on IZ and 7*, -;— 7*. Continuing this differentiation further at 

dt2 

each step i we shall obtain an equation of the form 

^(arHt2)7*(t2)) = ^{a^\t2h{t2))+j2UH,), 

where fij is a linear function of Hj with coefficients, depending on TZ 
and first i — 1 derivatives of j^, (this can be immediately seen by the 
induction). But at some stage the moments start to repeat themselves 
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due to equation (15. ip . So, taking K derivatives we shall obtain equa- 
tions of the following form 
(5.10) 

f (yi\t2)l,{t2) = a^\t2h{t2) + foo{Ho{t2)) 

j=o 



d^ ^ 



(it 2 (it 2 j=o 

-T7K{<^l\t2h.{t2)) = ^K^(t2)7(t2)) + E MHj) 

j=o 

Suppose that that the dimension of the inner state is n (i.e. dim?/ = 
n) , then we get that each of the matrices Hj has entries and there is 
the total number of n^(iV + 1) entries for the moments Hq, . . . , Hj^f. So, 
taking "enough" derivatives of 7*(t2) (i-e., taking K so that K dim{£y > 
n^{N + 1)) and eliminating all the entries of the moments, we shall ob- 
tain a finite number of polynomial differential equation for the entries 

of 7,(^2). □ 

Remark: From this theorem is follows that 7* satisfies an equation of 
the form 

± \ tAj J ilj J tij J , , , J JU \J ^ 

where P(xo, xi, X2, . . . , xk) is a non-commutative polynomial with co- 
efficients in 7?.. 



5.2. Sturm Liouville vessel parameters. The following example 
was extensively studied in |M2] . It deals with the Sturm Liouville 
differential equation 

d^ 

^2/(^2) - q{t2)y{t2) = >^y{t2), 

with the spectral parameter A. The parameter q{t2) is usually called 
the potential. For q{t2) = this problem is easily solved by exponents 
and in this case we shall call this equation trivial. In |M2] one connects 
solutions of the more general problem with non trivial 9(^2) to the 
trivial one. 
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Definition 5.7. the Sturm Liouville vessel parameters are given by 



di = 



1 

1 



(72 



1 




>7 



-tnnit^) -/3(t2) 

/3it2) I 





i 



for real valued continuous functions 7rii(t2), /?(^2)- 

The input compatibility differential equation (12. 161) is equivalent to 



a 

dt2 



Ui{X,t2) + iU2{X,t2 







Ml(A,t2) 
U2{X,t2) 



. From the second equation 



where we denote ma(^2) = 

one finds that M2(A,t2) = ^2) and plugging it back to the 

first equation, we shall obtain the trivial Sturm Liouville differential 
equation with the spectral parameter iX for ui{X,t2)- 

^Wl(A,t2) = lXUi{X,t2). 



For the output y\{t2) = 
equivalent to the system o: 



yi{X,t2) 

?/2(A,t2) _ 

equations 



we shall obtain that fl2.17p is 



(A - Z7rn(t2))2/i(A, t2) - + /?(t2))?/2(A, ^2) 
{f3{t2)-i-)yi{X,t2) + iy2{X,t2) = 



from which we immediately obtain that ?/2(A, ^2) = ^(/?(^2) — ^2) 
and plugging it into the first equation 



dt': 



;yi{X,t2) - {7fn{t2) + I3'{t2) + P\t2))yi{X,t2) = iXyi{X,t2 



which means that yi{X,t2) satisfies the Sturm Liouville differential 
equation with the spectral parameter iX and the potential ^'(^2) = 
(vrii(t2) + /3'(t2) + /3'(t2)). 



The first equation (11.21) considered for Hq 



-Z7rn(t2) -P{t2) 

from where we conclude that 
(5.11) H',' = -(3{t2), 



-"0 



ttU 7t12 
71721 7:t22 

-"0 -"0 



becomes 



TTii 77-22 
^0 - -"0 



-«7rii(t2). 
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Let US consider the differential equation (15 .4^ . where we use the nota- 



tion Hi — ^22 
parameters, we shall obtain 



Substituting the expressions for the vessel 



d TTl2 _ aul2 



0, 

Hl^-Hl\ 
Hp. 



and consequently, using the formulas (15. lip the second equation results 
in 



(5.12) vrn = — 

at2 

Together the first and the fourth equations give 



rii 
'0 



from where we obtain using f lS.lip 
(5.13) -^\m' + Hfi'] = Q 



dt2 



'0 



iTTiiH^' + 13 H^' + iH, 



t21 



Hi' 



Ce 



Additionally, Hq has to satisfy Hq 
shall obtain that 



^HqCi. Using this relation we 



Hi' 
Hi' 





"01" 






1 





m' 

HI' 



{HfY 

(Hir 
(H^r 
(Hir 



[Hi'r 



{Hfr 
(Hrr 





" 


1 " 




1 






from where we conclude that 



(5.14) 



H 



11 



iH 



22\* 
) 1 



H 



12 




:h 



12\* 




H 



21 



iH 



21\* 
) ■ 



As we can see H^^ 



/iq^ is a real valued (arbitrary at this stage) 



function, and Hr, is as follows 



(5.15) 



H^ 



r — ZTTii 



L21 

/4g 



-/3 
r + iTCii 



where r G M. 
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Let us perform the same calculations for Hi. From the algebraic equa- 
tion ([53]) 

we obtain that 

Hl^ + {Hiy Hl^ + {HfY 

From the equation (15. 4p with 2 = we obtain as before that 
(5.16) Hi' - H^ = ^Hl' + zvrni/o" + /^^oS 



-HqGi 'Hq 



dU 



(5.17) 



12 



d 

dto 



Hl'^ + iT^iiH^' + pH^' + iHl 



rl2 



22 



and the same equation (15. 4p with i = 1 produces similarly to the pre- 
vious case 

d 

dt2 

rll , tj22N 



(5.1^ 



d 



Hr-/3Hr + 2{Hl'-H!') = 0, 
-tT^iiHl' + ml'-HT)- 



(5.19) ^^{HV^Hi 

Plugging flHl^ and flETTD into <KYE^ . we shall obtain that H^' = h^' 
have to satisfy the following differential equation of the first order: 

d 



+tA.Hl'-T,,,Hl' + tPHl'=Q. 



or after cancellations 



(5.20) 2i{H^ 



2l\i 



2(7rni/o ) + +/3^ii^o " ^ii(^o " 



rll^ 



d? 



H] 



22 



dtr'' 



Inserting here the expressions for Hq appearing in ( I5.15P we shall ob- 
tain that the real part of the last equality can be derived from the 
equation ( I5.12p : 



0. 



The imaginary part gives the following equation 

r' 

Suppose (see |M2] ) that there exists a function r such that (3 = . 



Then using f l5.12p tth 



and inserting these equations into the 
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formula for -^h"^ we obtain that 



T 



{IV) 



T T 

Let us write down the formulas for the elements of Hi: 



t22 



22\ 



Hi ~ Hi - 



Hi 



dti 



dEUD 



The last equation f l5.20p ^ is obtained from (15.201) by substituting the 
index at Hq by the index 1: H^. 

Finally, we obtain that in the general case i/j+i is derived from Hi 
using a system of similar equations. Denote 



M. = i-Hf + ^nllHP + (5Hf + ^H^\ 
U = i-H^^+t7:iiHl' + (5Hf^ 



then 
(5.21) 



H. 



22 



TTl2 
TTll _ 

d ( TTll I 7t22 \ 

2^^^Si 



-Y + ^(^iiM)' 



from where we see that Markov moments are defined up to initial con- 
ditions for no = 2 elements. Notice that = 1 in this case. 
Let us also demonstrate theorem l5.6l for the Sturm Liouville parameters 
from definition 15. 71 We will take the simplest case n = 1 and as a result 
the transfer function is 



S{X,h 



1 



X + z 



C{h)B{t2)ai 



i=0 



-zyC{t2)B{t2)ai. 



We have already seen (in (15.121) ) that 7* is necessarily of the form 



7* (^2 



for a real valued function /3(t2) on I. In this case, = which means 
that the first moment is a multiple of the zero moment: Hi = —zHq. 
Since the vessel parameters are constant the differential ring IZ = C 
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is trivial. Using the formulas developed in |M2] we shall obtain that 
r = exp j 13 satisfies 

for k = \/—iz, which may obviously be rewritten as a polynomial 
differential equation for /3, after inserting the formula for r = exp J /3 
and multiplying by r^^. 

5.3. Partial realization problem. In this section we want to solve 
the partial realization problem 11.41 The key idea here is analyzing of 
the equations (11. 2p and fl5.4p . First, from the equation (11.21) we find 
that 

Consider now the system of equations (15.41) for i = 1, . . . ,n — 1, where 
we substitute 7*(t2) by the formula (11.21) above: 

- a^^j^Hi + Hia^^-f = erf V2i^i+i - i^i+iaf V2 
= erf V2i/i+i - i/i+i(TfV2 

= --^Hi + {a2Ho - aiHQa^^a2)Hi + af V2i/i+i - i/i+ifff V2 
This condition can be rewritten in a more compact way as 

(5.22) [Hi,a^'j] = -^H, - ai[Ho, a^'a2]H, - [H,^i, a^'a^], 

at2 

where [A, B] = AB — BA denotes the commutator of A and B. 
So, the question which arises here is whether every function 7, which 
solves (I5.22P produces a solution of the partial realization problem 11.41 
We connect the solution of the partial realization problem 11.41 to the 
classical one. 

Theorem 5.8. Assume that we are given Junctions Ox^a^ and Markov 
moments Hi(t2),i = 0, . . . ,n (n > 0) defined in a neighborhood of a 
point ^2 £ I; satisfying the necessary conditions (15. 5 p and (15. 6p . As- 
sume also that the following necessary conditions are also satisfied: 

(1) The classical partial realization problem for the moments Hi{t^ 
i = 1, . . . ,n has a solution, which we denote by S & TZS with a 
minimal realization 

S{X, tl) = I- BoX^\XI - Ao)-'Bo(Ji{tl) 
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(2) There exists a solution 7(t2) of the equations f l5.22p for each 
i G {0,...,T2,}, and such that B{t2), X(t2) defined from the 
formulas (14. ip . (14. 3p using Bq, Aq, Xq and 7(^2) realize the 
first moment: 

Ho{t2) = {B{t2)yX{t2)-'B{t2)ai{t2). 

Then, the function 

S{X, t2) = I- Bit2)X{t2)-\XI - Ao)-^B{t2)a,{t2) 

is a solution of the partial realization problem \1.4\ with 7*(t2) given by 
the equation (11.21) . 

Proof: Using simple calculation and the definitions of i?(t2),X(t2) it 
follows that the function S{X, ^2) intertwines solutions of the input and 
output LDEs, with 7^, defined as in the theorem, is equal to Ip a.t 00, and 
is continuous for each A. Since the Lyapunov equation holds, the func- 
tion ^(A, t2) G TZSX. It remains to show that it has the given Markov 
moments. By construction the first moment of S{X,t2) is ifo(^2)- By 
Theorem 15.31 the second moment of S'(A,t2) satisfies (15. 4p with i = 0, 
and differential equations with hq initial conditions. But these initial 
conditions are obtained from the realization S{X,t2). This way we ob- 
tain the first moment of ^(A, ^2)- The other moments are obtained by 
iteration. □ 

Remarks: 1. If there is a linear combination -^(^2) = X] cii-^j(^2) 

which has its spectrum disjoint from —H{t2) then 7 is uniquely deter- 
mined from 

n-l 

^ ra— 1 n— 1 n— 1 

= --7- Yl OiiHi - ai[HQ,a^^a2\ Yl (^i^i - [Yl Oi-f^i+i, o-f V2] 

2. Case cri(t2) > 0: The equation (I5.22p is then uniquely solvable 
for 7 if the spectrum because the first Markov moment .^^0(^2) is self 
adjoint and strictly positive, and thus the spectra of Hq and —Hq are 
disjoint. 

6. Nevanlinna-Pick interpolation problem 

Nevanlinna-Pick (NP) interpolation problem 11.61 is similar to non com- 
mutative, Riemann surface cases and involves specifying the exact class 
of input-output mappings, and a finite number of inputs that are to 
be mapped to the corresponding outputs. In our case this yields fixing 
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both U(A, t2) and U*(A, — — ; t2), a set of values Wj, i = 1, . . . , for 

the spectral parameter, and corresponding solutions ^4(^2) and 774(^2) 
of the input and the output LDEs, respectively, as it is defined in the 
problem 11.61 

Since everything is determined by the initial conditions, this reduces to 
the NP problem at some fixed G I, except that now 7* is fixed, yield- 
ing a constraint on the solution parameter Sq. According to remark, 
following the theorem 15.61 7* satisfies an equation of the form 

for a non-commutative polynomial with coefficients in IZ. Since the 
derivatives of 7^, may be represented using the moments Hi{t2) differ- 
entiating the linkage condition fl2.8p and using f l5.4p . we shall obtain 
that the moments have to satisfy an equation of the form 

P'(i/o,i/i,...,i/„) = 

for a non-commutative polynomial P' with coefficients in 7?.. On the 
other hand, they must be linearly dependent, which is a necessary con- 
dition and the formula (15.11) must hold. Consequently, there must be a 
minimal (with respect to the index - the biggest moment appearing) 
linear combination 

N 
i=l 

of the moments such that it divides from the left and from the right 
the polynomial P'{Hq, Hi, . . . , Hn). In this way we shall be obtain the 
coefficients /ij, which are also coefficients of the minimal polynomial for 
the final operator Ai. 

Next we shall look for functions Sq G IZS, which has Ai with the 
prescribed minimal polynomial. Notice that this also means that the 
Jordan block structure of the class of Ai is defined in this manner. It 
is remained to choose values for Bq, Xq so that the system of equations 
(I5.10p will be satisfied at the point ^3- This follows from the fact that 
7*^0 corresponding to the function 5*0 will have the same initial values as 
the given 7* and, moreover, there will be obtained the same differential 
equations for the 7=,, and 7*^0- We now present the exact statement 

Theorem 6.1. Suppose that there is a realization of Sq G IZS 

So{X,t2) = I - B*X-\XI - Ai)-iSoai(t°), 

then So is a solution of the Nevanlinna-Pick interpolation problem \1.6[ 
if the minimal polynomial of Ai is equal to the minimal polynomial 
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defined from 7* and if the moments Hi of the functions Sq satisfy the 
system of equations (15.1 Op at 

Proof: We construct first a polynomial P{x,x',x", . . . ,x^^^) from 7* 
as it is explained in the text preceding the theorem. Then we find a 
minimal polynomial for Ai, obtained from P'[Hq, Hi, ... , H^) by sub- 
stituting Hi with X*. Then the moments of the function S'o(A, t2)'K-SX 
constructed from Sq G 7?.«S will satisfy the same differential equations 
appearing in theorem flS.lOp . As a result the output 7*^0 of the function 
So{\it2) will satisfy the same polynomial differential equation as the 
given 7* and will have the same initial conditions, which means that 
7*,o = 7*) which finishes the proof. □ 
We now solve Problem II. 7[ 

Theorem 6.2. Given C^^^-valued functions cji, 02, 7, an interval I 
and N quadruples < 4, Wj, ^j, rjj > where ^2 ^ I, Wj G C_|_, ^j, rjj G C^^^ 
j = 1, . . . N , and assume that the corresponding matrices Xj > 0. Then 
there exists a solution of the Nevanlinna-Pick problem i.e. there 
exists a function S G IZSX satisfying S{wi,t2)Ci = Vi ^'^^ only 
if there exists n G N, matrices A*),X^,G C^"-^)''^""^) with X^ > 0, 
g Cpy-in-i) ^y,, ^ (^nxn ^^^/^ ^/^^^ ^.^ q,^ ^, ^^fined by 



(6.1) 
(6.2) 



Bi 







Vi - ii 



x^ 





X,- 



(6.3) Ai 



-7^ia,{t\){BlY{^^y 



X, 



X,: 



''2 5 



4)dA 



VijBj. 



it holds that 

(1) A, = V,,A,V,J\ 

(2) §{XI-A,)-^Bia^\ti)<l>-\X,t 
and the matrix X{t2), 

X{t2) = X, + r B,{y)a2{ym{y)rdy 

is invertible on the interval I. 

Proof: For each ^2 all the functions which satisfy S{wi,t2)^i — Vi 
of the form Te-(5^(A, ta)), provided X, = ^^^^('^2)4* - ^.^1(^2)^; 

W* + Wi 



> 0. 
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Here 



e.(A) 



/ + 



X(A 



and Sn G "K-S. Given a minimal realization 



-\XI-A^,r'Bla,{i^,) 



of 5*0, we shall obtain from formulas f l3.15p . f l3.16p . f l3.17p the formulas 
fl6.ip , fl6.2l) , fl6.3p in the theorem. In view of Theorem 14.71 a necessary 
and sufficient condition to obtain the same function S{X,t2) for every 
i is that there exist invertible constant matrices Vij such that the op- 
erators Ai are similar. In other words there must exist Vij such that 
Ai = VijAjV^^. Moreover, the second part of theorem 14.71 tells that 
additionally the equality 



{XI - A,)-'B,a^\ti)<l>-\X,tif^)dX = VjBj 



must hold. 



□ 
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ON THE CLASS SX OF J-CONTRACTIVE 
FUNCTIONS INTERTWINING SOLUTIONS OF 
LINEAR DIFFERENTIAL EQUATIONS 



DANIEL ALPAY, ANDREY MELNIKOV, AND VICTOR VINNIKOV 

Abstract. In the PhD thesis of the second author under the 
supervision of the third author was defined the class SX of J- 
contractive functions, depending on a parameter and arising as 
transfer functions of overdetermined conservative 2D systems in- 
variant in one direction. In this paper we extend and solve in the 
class SX, a number of problems originally set for the class S of 
functions contractive in the open right-half plane, and unitary on 
the imaginary line with respect to some preassigned signature ma- 
trix J. The problems we consider include the Schur algorithm, 
the partial realization problem and the Nevanlinna-Pick interpo- 
lation problem. The arguments rely on a correspondence between 
elements in a given subclass of SX and elements in S. Another 
important tool in the arguments is a new result pertaining to the 
classical tangential Schur algorithm. 
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1. Introduction 

Functions S{X), which are C''^^^- valued, analytic and contractive in the 
open right half plane C_|_, or equivalently, such that the kernel 

X + w* 

is positiv^ in C+, play an important role in system theory, inverse 
scattering theory, network theory and related topics; see for instance 
[L] . |BC] . |DD] . [He], Here, positivity of the kernel means that for 
every n G N and every choice of points Wi, . . . ,Wn G C+ and vectors 
• • • ) £ C^^^ the n X n Hermitian matrix 

is positive (that is, has all its eigenvalues greater or equal to 0). 

Far reaching generalizations of this class were introduced in [Ml iMVTl 
IM Vc] ■ in the study of 2D-linear systems (say, with respect to the vari- 
ables (^1,^2)), invariant with respect to the variable ti. To introduce 
the classes defined in these papers we first need a definition. 

Definition 1.1. Let ai, 02, 7 and^^, be C^^''^ -valued functions, contin- 
uous on an interval I = [a,b]. Suppose moreover that a\ and take 
self-adjoint values, and that ai is differentiahle and invertible on I, and 
that the following relations hold: 

7(^2) + 7(^2)* = 7*(^2) + 7*(i2)* = -^(Tiit2), t2 e I. 

dt2 



^positive-definite in the classical terminology 
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Then cr2, 7, 7=,, and the interval I are called vessel parameters. 

The class of functions SI corresponding to thr vessel parameters ai, 
(T2, 7, 7* and the interval I was introduced in jM | IMVc] (see Definition 
12.71 below) and consists of the functions S{X,t2) of two variables X,t2 
such that for every ^2 ^ I the function S{X,t2) is meromorphic in C+ 
and the kernel 

^ ' ^ X + w* 

is positive for A and u in the domain of analyticity of S{X,t2) in C4.. 
For positive o'i{t2), the positivity of the kernel implies that S is analytic 
in C+; see |Do] . [X] . For general (invertible) (Ji{t2), the entries of 5* are 
of bounded type and S has at most poles in C+; see |ADRSj . It is also 
required that S{X,t2) is analytic at infinity for each t2, with value Ip 
there, and that S{X,t2) maps solutions of the input Linear Differential 
Equation (LDE) with the spectral parameter A 

d 

Xa2it2HX,t2) - aiit2)—u{X,t2) + -f{t2)uiX,t2) = 0, 

Ot2 

to solutions of the output LDE 

d 

Xa2{t2)y{XM) - (yi{t2)Tz-y{XM) + l*{t2)y{X,t2) = 0. 

012 

It is proved in [Ml IMVcj that elements of SX are the transfer functions 
of ti-invariant conservative 2D systems; see Section [21 

The purpose of this paper is to study various questions in the class SI 
such as Nevanlinna-Pick interpolation, property of moments, etc. A 
key result is the following theorem, which we prove in the sequel; see 
Section HI 

Theorem 14.11 Let us fix the parameters ai,o'2, and 7, and the interval 
I. Then for every t2 & I there is a one-to-one correspondence between 
pairs (7*,5') such that S G SX and continuous in a neighborhood 
of t\, and functions Y{X), meromorphic in C+, and with the following 
properties 

(1) r(oo) = Ip, 

(2) Y{X)*(yi{tl)Y{X) < aiitl) for A G C+ where Y zs analytic, and 

(3) r(A)Vi(t^)r(A) = CTi(t^) for almost X satisfying ^X = and 
where Y{X) is the non-tangential limit. 

As mentioned above the (Ji(t2)-contractivity of Y implies that Y is of 
bounded type in C+, and thus the asserted non-tangential limits exist 
almost everywhere. 
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Definition 1.2. The class of functions Y with the properties in Theo- 
rem \4.1\ will be will be denoted by Slt^) and the functions will be called 
ai{t^-inner. The subclass of rational functions of S{t'^ will be denoted 
by TiS{tl). 

For the sequel, it is important to notice that the general tangential 
Schur algorithm developped in |AD] can be applied to functions in 
«S(t2), and in particular in T^SiJ^l)- 

Definition 1.3. IZSX will denote the subclass of functions in SI 
which are rational in A for every t2 

The paper consists of six sections besides the introduction, and we 
now describe their content. In Section [2] we review the main results 
from [Mj (see details in |MVlj and |M Vcj ) on ti-invariant conservative 
2D-systems, relevant to the present work. In particular the class SX 
mentioned above consists of the transfer functions of these systems. 
In Section [3] we present the reproducing kernel space approach to the 
tangential Schur algorithm for the class <S, as developed in |AD] . We 
obtain in particular new formulas which allow us to find the main op- 
erator in a realization of an element of «S after one iteration of the 
tangential Schur algorithm; see formulas fl3.15p . fl3.17p . f l3.16p in The- 
orem 13.51 In Section S] we develop the tangential Schur algorithm for 
a function S{X,t2) G SX. Applying directly the theory of the previous 
section to 5'(A, ^2) leads to a new function which need not belong to 
SX. Instead, we apply the Schur algorithm to the (Ti(t2)-inner function 
S{X, tg) for some preassigned ^2 ^ I? and obtain a simpler (in terms of 
McMillan degree) ai(t2)-inner function So{X,t2). We use Theorem 14.11 
to obtain an element in a class SX from So{X,t2). We call this pro- 
cedure the tangential Schur algorithm for the class SX. We study in 
Section |5] the coefficients (called the Markov moments) Hi{t2) of the 
expansion of S{X,t2) around A = 00 

00 ^ 

(1.2) s{XM) = Iv-Y.^lH^ii2) 

i=0 

It turns out that the first Markov moment //o(^2) satisfies the Lyapunov 

equation 

(1.3) 

7*(^2) - 7(^2) = (T2{t2)Ho{t2) - a,{t2)Ho{t2)a{\t2)a2{t2), t2 e I, 

which means that given the functions cri,cr2,7 and Hq on I, one can 
uniquely reconstruct 7,,, and, as a result of Theorem 14.11 there will 
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exists a unique function S'(A, ^2) {^2 £ I) with the given Markov param- 
eters. We prove the following theorem 15.51 on the structure of Markov 
moments: 

Theorem Suppose that we are given moments Hi{t2), defined in a 
neighborhood of the point ^2 ^ I- Then there exists Uq < such that 
each element i^j+i is uniquely determined from HQ,...,Hi using the 
algebraic formulas (15. 5p . andriQ LDEs with arbitrary initial conditions, 
obtained from (15.41) . Moreover, Uq equations must be satisfied by the 
elements o/7*(t2)- 

Starting from the Markov coefficients at infinity, constructed under 
conditions of Theorem 15.51 we recover in Section E] the transfer function 
S{X,t2). In the arguments we make use of a Krein-space realization 
theorem of Dijksma-Langer-de Snoo for analytic functions at infinity 
[DLdeSl Theorem 3]. By a counterexample, we show that it is not 
always possible to reconstruct the function using Hilbert space realiza- 
tions. 

In Section [7] we study two generalizations of Nevanlinna Pick inter- 
polation problem in the class SX. Let us recall that in the classical 
Nevanlinna Pick interpolation problem [Nl |P] Schur analysis plays a 
special role; see for instance [FKj . |Dy| , [X] . Schur analysis gives a 
parametrization of all solutions (when they exist) for the given data. 
The first problem 11.41 is similar to the non commutative and Riemann 
surface cases and involves specifying the exact class of input-output 
mappings, and a finite number of inputs that are to be mapped to the 
corresponding outputs. In our case this yields fixing vessels parameters 
and a set of values Wi,i = 1, . . . , N and corresponding solutions Ciih) 
and r]i(t2) of the input and the output LDEs, respectively: 

Problem 1.4 (Nevanlinna-Pick interpolation) . Lei cti, (J2, 7, 7* be ves- 
sel parameters and let 1 be an interval. Let N eN and Wi,i = 1, . . . , N 
be complex numbers. Suppose also that N input functions ^i{t2) satis- 
fying (12.161) with corresponding spectral parameters Wi 's, and N output 
functions r]i{t2) satisfying (12.171) with the given Wi's. 

(1) Give sufficient and necessary conditions, so that there exists 
S{X,t2) e SI such that S{wj,t2)^i{t2) = Viih), j = l,---,^; 
on a sub-interval of I. 

(2) Describe the set of all solutions for this problem. 

The second problem uses the fact that we can also specify the data for 
different values of t2 and is more similar to the classical one: 
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Problem 1.5. Given -valued functions cri,(T2,7 defined on \, and 
given N quadruples < tljWjy^jjTjj >, where 4 £ ""^i ^ ^j^Vj ^ 
C^^P J = 1,...N, then: 

(1) Give sufficient and necessary conditions, so that there exists 7* 
and S{X, ^2) £ <5X such that S{wj, = Vj> J = 1) • • • ; '^'^ 
a sub-interval of I containing all the 4 ■ 

(2) Describe the set of all solutions for this problem. 

If all the values t^ = ^ are equal, we have to find a function S{X, t^) 
satisfying S{wj, t^)^- = Vj- Thus, the above problem is a generalization 
of the classical Nevanlinna-Pick interpolation problem. We also remark 
that we do not address the question of describing the set of all solutions. 
Remarks: Some of the results presented here have been announced in 



2. ti INVARIANT CONSERVATIVE 2D SYSTEMS. 

The material in this section is taken from jM| IM Vll IM Vc] . where proofs 
and more details can be found. The origin of this theory can be found 
in the paper [Lij. 



2.1. Definition. An overdetermined conservative ti-invariant 2D sys- 
tem is a linear input-state-output (i/s/o) system, which consists of 
operators depending only on the variable ^2 and is of the following 
form: 



(2.1) /S : < 



J^x(ti,t2) = Ai{t2)x{ti,t2) + Bi{t2Mti,t2) 

x{ti,t2) = F(t2,t2)a;(^i,4) + / F{t2,s)B2{s)u{ti,s)ds 



[ y{ti, t2) = uitiM) - B{t2yx{ti,t2), 

where the variable ti belongs to M, and the variable ^2 belongs to some 
interval I. Furthermore, the input u{ti,t2) and the output y(ti,t2) 
take values in some Hilbert space S and the state x{ti,t2) takes values 
in the Hilbert space T-Lt2- We assume that u{ti,t2) and 2/(^1,^2) are 
continuous functions of each variable when the other variable is fixed. 
The operators of the system are supposed to satisfy the following: 

Assumptions 2.1. 

(1) Ai{t2) : Ht^ — ^ 'Ht2, and -8(^2) : £ "Htj o-re bounded operators 
for all t2, 

(2) The functions di, o"2, 7, 7* : — )■ are continuous in the oper- 
ator norm topology. 
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(3) 0"i(t2) is an invertible operator for every t2 G I. 

(4) F{t, s) is an evolution continuous semi-group. 

For continuous inputs u(ti,t2), the inner state is continuously differ- 
entiable. Requiring now the invariance of the system transition from 
(ti, to (ti, t2) via the points (t?, ^2) and (ti, ^2) respectively, is equiv- 
alent to the equality of second order partial derivatives of ^2): 

Substituting in this equality the system equations we obtain that for 
the free evolution u{ti, ^2) = the so called Lax equation holds 

(2.3) Ar{t2) = F{t2A)Mtl)F{tlt2)- 

Inserting fl2.3p into fl2.2p we see that the input u{ti,t2) has to satisfy 
the following PDE 

B{t2)cT2{h)i-^u{ti,t2) - B{t2)ai{t2)i-^u{ti,t2)- 

{A,{t2)B{t2)a2it2) + F{t2,tl)^^[F{tlt2)B{t2)a,{t2)])u{t,,t2) = 0. 

Assuming the existence of a function 7(t2) satisfying 
(2.4) 

Al{t2)B{t2)cr2{t2) + F{t2,tl)-^[F{tlt2)B{t2)ai{t2)] = -B{t2h{t2) 

OS 

we obtain that it is enough that u{ti, ^2) satisfies the PDE 

(2.5) a2{t2)^u{ti, t2) - (yi{t2)^u{tiM) + 7(^2)^(^1, ^2) = 0. 

The output y{ti,t2) should satisfy the output compatibility condition of 
the same type as for the input compatibility condition (12. 5p . namely: 

(2.6) a2(t2)^2/(ti, ^2) - cri{t2)—y{t,,t2) + 7*(^2)2/(ti, ^2) = 0. 

Inserting here y(ti,t2) = u(ti,t2) — B{t2)*x{t2,t2) we obtain that 

(2.7) = a2{t2)B{t2yA,{t2)F{t2,tl)- 

-ai(t2)T^[5(t2)*F(t2,t^)]+%(t2)S(t2)*F(t2,t^) 

012 

(2.8) 7(^2) = criit2)B{t2yBit2)cj2it2)~ 

~(J2{t2)B{t2yB{t2)ai{t2)+l.{t2). 
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The fact that the system is lossless comes from the requirement of the 
so called energy balance equations: 

d 

— (x(ti,t2),a;(ti,t2))%2 + {(^My{ti,i2),y{ti,t2))£ = 

= {cri{t2)u{ti,t2),u{ti,t2))£, i = 1,2, 

which means that the energy of the output is distributed between the 
energy of the input and the change of the energy of the state of the 
system. Immediate consequences of this requirement are 

(2.9) = Aiit2) + Alit2) + B{t2)a,it2)B{t2)*, 

(2.10) -^[F*(t2, f2)F{t2, i)] = F*{t2, 4)B{t2ya2it2)B{t2)F{t2, 4)- 
at2 

In this manner we obtain the notion of conservative vessel in the integral 
form, which is a collection of operators and spaces 

QJ = (A(t2), F{t2, B{t2); ai{t2), ^2(^2), 7(^2), 7*(^2); nt„S) 



where the operators satisfy the regularity assumptions 12.11 and the 
following vessel conditions: 

= A,{t2) + Al{t2) + B{t2yMt2)B{t2) dH 

||F(t2,t^Mti,t°)|p-|Kti,tO)f =^ 

= j',f{a2{s)B{s)x{ti,s),B{s)x{ti,s))ds I^M) 

F(t2,t^)Ai(tO) = Ai(t2)F(t2,t^) m 

= ^(F(tO,t2)S(t2)ai(t2)) + 

t2)Ai{t2)B{t2)a2{t2) + F{tl ^2)5(^2)7(^2) dH 

= ai(t2)^[^(t2)*/(t2,t^)]- 

-(y2{t2)B{t2yA^{t2)F{t2A)-l*{t2)B{t2rF{t2A) (D 

7(^2) = -02{t2)B{t2yB{t2)ai{t2) + 

+(Tl{t2)B{t2yB{t2)(T2{t2) + 7*(^2) (ES]) 

In order to simplify some notations we introduce the following defini- 
tion. 

Definition 2.2. Let 

U = U(Ai, A2; t2) = a2(t2)Ai - ai(t2)A2 + 7(^2), 
and similarly 

U* = U*(Ai, A2; ^2) = o-2(t2)Ai - o-i(t2)A2 + 7*(^2)- 



ON THE CLASS SI 

The vessel QJ is naturally associated to the system (12. ip 

dti 



S : i 



t2 _ 

iti,t2) = F{t2,t'^2Mti,4) + I F{t2,s)B{s)a2{s)u{ti,s)ds 



^ y{tiM) = u{h,t2) - B{t2Y X{tut2). 

with inputs and outputs satisfying the compatibility conditions f l2.5p 
and fl2.6p . i.e. satisfy: 

U(7^' h>{tut2) = 0, U,( — , — ; t2)y{h, ta) = 

Oti Ot2 Oti Ot2 

The theory of such vessels, developed in [Ml IMVl] enables to find a 
more convenient form of the vessel. Denoting = "Hto, A\ = Ai{t2), 

F*{t2,4)F{t2,tl) = X-\t2) andS(t2) = ^2)5(^2), we shall obtain 
the following notion, first introduced in |M2] . 

Definition 2.3. A (differential) conservative vessel associated to the 
vessel parameters is a collection of operators and spaces 

(2.11) QJ= (Ai,S(t2),X(t2);ai(t2),a2(t2),7(^2),7*(^2);'H,£), 
where the operators satisfy the following vessel conditions: 

(2.12) = A(5(t2)ai(t2)) + A,B{t2)a2it2) + B{t2)i{t2), 

dt2 

(2.13) A,X{t2) + X(t2)At = B{t2)a^{t2)B{t2)\ 

(2.14) 4-X(t2) = 5(^2)^2(^2)5(^2)*, 

dt2 

(2.15) 7*(i2) = 7(^2) + cy2{t2)B{t2y^~\t2)B{t2)(y^{t2)- 

-Mt2)B{t2yX-\t2)B{t2)a2{t2) 

It turns out that Lyapunov equation (12.130 is partially redundant. 

Lemma 2.4. Suppose that B{t2) satisfies (I2.12p and li{t2) satisfies 
(I2.14p . then if the Lyapunov equation (I2.13P 

A^X{t2) + X(t2)At + Bit2)(x^it2)B{t2y = 

holds for a fixed tg? then it holds for all ^2- U ^{t^) — X*(t2); then 
X(t2) = X(t2)* for all t2. 

Proof: By differentiating the left hand side of (I2.13p . we will obtain 
that it is zero. The derivative of X(t2) is selfadjoint. □ 
This representation of a vessel is the most convenient when one focuses 
on the notion of transfer function, as we do in the next subsection. 
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2.2. Transfer function. Performing a separation of variables as fol- 
lows 

U{ti,t2) =ux{t2)e^'\ 

y{tiM) = yx{t2)e^'\ 

we arrive at the notion of a transfer function. Note that u{ti,t2) and 
y{ti, ^2) satisfy PDEs, but U\{t2) and yx{h) are solutions of LDEs with 
spectral parameter A, 

(2.16) \J{X,^;t2)u{h,t2) = 0, 

012 

(2.17) U(A,^;t2)l/(ti,t2) = 0. 

The corresponding i/s/o system becomes 

f xxit2) = (A/ - A,{t2))-'B{t2)aiit2)uxit2) 

i-xx{t2) = F{t2,4)xxitl) + J F{t2,s)B2{s)ux{s)ds 

^ yx{t2) = Ux{t2) - B{t2yXx{t2) 

The output yxih) = uxih) — B{t2)*xx{t2) may be found from the first 
i/s/o equation: 

y\{t2) = S(\,t2)ux(t2), 
using the transfer function 

S{X,t2) =I-B{t2nXI~A,{t2))-'B{t2)ai{t2) 

= I - S(t2)*(A/ - F{t2, t°) A(tO)F(tO, t2))-'B{t2)a,{t2) 
= I - B{t2yF{t2,tl){XI - A^{tl)y^F{tlt2)B{t2)a^{t2) 

= I - B{t2yF*{tlt2)F*{t2A)F{t2A){^I - A,)-^B{t2)cy,{t2) 
= I - B{t2yX-\t2){XI - A,)-^B{t2)a,{t2) 

and we obtain that 

(2.18) 5(A, t2) = I- B{t2yX-\t2){XI - A^)-'B{t2)a^{t2). 

Proposition 2.5. The transfer function S{X,t2) defined by (12.181) has 
the following properties: 

(1) For allt2, S{X,t2) is an analytic function of X in the neighbor- 
hood of 00, where it satisfies: 

S{oo,t2) = Ip 

(2) For all X, S{X,t2) is a continuous function 0/^2- 
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(3) For A in the domain of analyticity of S{X,t2) : 

(2.19) Six, ta) Vi(t2)^(A, h) < a^ih), > 0, 
and 

(2.20) S{\ t2yai{t2)S{\ t2) = Mt2), 3?A = 0. 

(4) Maps solutions of the input LDE fl2.16p with spectral parameter 
A to the output LDE (12.171) with the same spectral parameter. 

Proof: These properties are easily checked, and follow from the defi- 
nition of S'(A, ^2): 

S{\ t2) = I- B{t2nXI - A,{t2))-'B{t2)a^{t2). 

The function S'(A,t2) is analytic for A > ||yli(t2)|| and since all the 
operators are bounded, we have S{oo,t2) = Ip. The second property 
follows from the regularity assumptions 12.11 The third property follows 
from straightforward calculations: 

5(A,t2)Vi(t2)5(A,t2)-ai(t2) = 

-2^{X)a^{t2)B{t2yCXI - Al{t2))-\XI - A^{t2))-'B{t2)cr^{t2) 

Here the sign of — 3?(A) determines the sign of '5'(A, t2)*c"i(^2)5'(A, ^2) — 
(71(^2) and thus the third property is obtained. The fourth property 
follows directly from our construction. □ 

Remark: When dim'H < 00, we obtain that S{X,t2) is a rational 
function of A for every t2- 

It is an interesting fact that also the converse of Proposition 12.51 holds; 
see [M], [MVcl chapter 5]. 

Theorem 2.6. For any function of two variables S{X,t2), satisfying 
the conditions of Proposition l275\ there is a conservative ti-invariant 
vessel whose transfer function is S{X,t2)- 

We define the class of transfer functions mentioned in the introduction 
as follows: 

Definition 2.7 f pVc] ). The class SI = 5X(U(A, — ; t2), U,(A, — ; ^2 

OT2 (yt2 

consists of functions S{X,t2) of two variables, which 

(1) are analytic in a neighborhood of X = 00 for all t2 and where it 
holds S{oo, ^2) = Ip, 

(2) are continuous as functions of ^2 for all X, 

(3) satisfy (12.191) and (12.201) in the domain of analyticity of S, 



12 



D. ALPAY, A. MELNIKOV, AND V. VINNIKOV 



(4) map solutions of the input LDE ^2.10i) with spectral parameter 
A to the output LDE l\2.17\ ) with the same spectral parameter 

Recall (see |CoLe] ) that to every LDE can be associated an invertible 
matrix (or operator) function $(^2, ^2)5 called the fundamental solution, 
which takes value / at some preassigned value cind such that any 
other solution u{t2) of the LDE, with initial condition "^(^2) = ^0 is of 
the form 

U{t2) = $(t2,t2)%- 

Let <l'(A,t2,^2) '^'*(A, ^2, ^2) t>e the fundamental solutions of the 
input LDE fl2.16p and the output LDE fl2.17p respectively, where we 
have added in the notation the dependence in A. Then, 

(2.21) Six, t2)$(A, t2, tl) = <I>,(A, t2, t°)>S(A, t°) 

and consequently 5'(A, ^2) satisfies the following LDE 

(0 00) i;SiX,t2) = a^\t2)ia2it2)X + -f.it2))SiX,t2)- 

^ ' -Six, t2)a^\t2)ia2it2)X + 7(^2)). 

The symmetry conditions f l2.19p and (12 .200 . which are a result of Lya- 
punov equation fl3.8p gives us 

(2.23) SiX,t2) = a^\t2)S-^*i-X,t2)a,it2) 

3. SCHUR ANALYSIS FOR THE CLASSICAL CASE 

One of the approaches to the tangential Schur algorithm for a rational 
matrix function S'(A) G C^^^ is based on the theory of reproducing ker- 
nel Hilbert spaces of the kind introduced by de Branges and Rovnyak; 
see [dBRlj . |dBR2] . Dy for information on these spaces. The paper 



[ADj considers the case of column-valued functions. In this section 
we adapt the results of [AD] to the case of row-valued functions. We 
also present a new formula for a realization of a Schur function after 
implementation of the tangential Schur algorithm. 

3.1. Schur functions and Reproducing Kernel Hilbert spaces. 

In this section ai denotes a fixed self-adjoint and invertible (but not 
necessarily unitary) matrix in C^^^. Let S'(A) be a rational function, 
cji-inner in the open right half plane, i.e. 

5(A)Vi5(A) - cTi < 

at all points in the domain of analyticity fi(S') of S in C+, and 

^(A) Vi,S(A) - ai =0 
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at all points on the imaginary axis where S is defined. Then, the kernel 

ai - S{wyaiS{X) 



{3.1] 



w + X 

is positive for X,w E ^{S), and the space of rational C^^^-valued func- 
tions 

n 

1=1 

is finite dimensional. These well-known facts can be proved using re- 
alization theory; see for instance [AGj . Furthermore, 'H(S') is the re- 
producing kernel Hilbert space, associated to the kernel -^^(A, w). The 
inner product is defined by 

{cKsiX,u),dKsiX,w))ns = {cKs{w,p))c^xp = cKs{w,iy)d*. 
For an arbitrary /(A) G 'H(S') we have the reproducing kernel property 

{f{X),^K{X,w))ns = f{w)C- 

More generally, let now be a finite dimensional Hilbert space of 
C^^^-valued functions defined in some set Q, and let {/i(A), . . . , fN^X)} 
be a basis of Ai. Let X G C^^p denote the Gram matrix with £,j entry 
given by 



(3.2) 



(/,(A),/KA))^, i,j = l,...,p. 



It is easily seen that the space is a reproducing kernel Hilbert space 
with kernel given by the formula 

, /i(A) 

(3.3) K{x,w)=[Mwy ■■■ Mwy]x-' 



I fN{X) 



We set 
(3.4) 



F(A) 



/i(A) 



L /p(A) 

Assume now that A4 consists of rational functions, defined on a set 
fl{^A). For a G the backward-shift operator Ra is defined by 

X — a 

Suppose that: 

(1) The space Ai is invariant under the action of Ra, 
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(2) The functions /j(A) has the property that /i(oo) = 0, i.e., 
F(oo) = 0, 

then the function F given by (13. 4p can be written as 

F(A) = {XI-A)-^Bai, 

for suitably chosen matrices A, B. In this special case, formula (13. 3 p 
takes the form 

(3.5) K{X, w) = aiB*{wI - A*)-^X-\XI - A)-^Bai. 

As mentioned at the beginning of this section we are interested in 
kernels of the form (13. ip for some ai-inner rational function S. We now 
recall the characterization of these spaces, and first note the following: 
equation (13. ip leads to 

(3.6) - SiwYa^SiX) ^ ^^^y^-,^^^y 

w + X 

If S is analytic at infinity and satisfies there S{oo) = Ip, and letting 
— )■ oo in this equation, we obtain the formula 

(3.7) S{X) = Ip- B*X-\XI - A)-^Bai. 

Theorem 3.1. Let A4 be a finite dimensional Hilbert space of V}^^- 
valued rational functions, which are zero at infinity. Suppose that 
Ra-M. C M. for a G Vt{M.) and let X he its Gram matrix with re- 
spect to F{X). Then M. = l-L^S) for S defined by (13. 7p if and only if 
the Lyapunov equation 

(3.8) AX + XA* + B(TiB* = 
holds. 

When the spectrum of the operator A is in the open left half plane 

3? A < 0, one has 

H(5) = H2,^,eH2,.,5, 

where ii2,cri is the Hardy space with the inner product 



We set J 



-ai 
ai 



Note that J is both invertible and self- 



adjoint, and one can define J-inner rational functions. Let G 



©21 ©22 

be a J-inner rational function; we introduce the linear fractional trans- 
formation 

Te{W) = (On + I^e2i)-'(ei2 + ^^622). 
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Theorem 3.2. Let S and be respectively ai-inner J -inner rational 
functions. Then there exists a ai-inner rational function W such that 
S = Tq{W) if and only if the map 



(3.9) 



F 



-S{X) 
L 



is a contraction from 'H(0) to 'H(5'). 

This theorem originates with the work of de Branges and Rovnyak (see 
|dBRlt Theorem 13 p. 305]), where it is proved in a more general 
setting, and is one of the key ingredients to the reproducing kernel 
approach to the Schur algorithm. The proof is the same as for column- 
valued functions, and is omitted. 

As a special case of the previous theorem we have: 

Corollary 3.3. Given S G IZS, Wi,...,Wn G C+, and row vectors 

B= ■■ ■■ , Ai = dmg[-wl,...,-wl]. 

Let M be a Hilbert space of row vectors spanned by the rows of the 
matrix-valued function 

F{X) = {XI-Ai)-^BJ. 
Let X be the solution of the Lyapunov equation 

(3.10) AX + XA* + BJB* = 
and assume X > 0. Let G be defined by 

(3.11) e(A) = hp - B*X-\XIn - Ai)-^BJ. 
Then there exists So E TtS such that S = Te(5'o). 

Proof: Let us denote by t the map ( 13. 9p . For an arbitrary element 
/(A) = r]F{X), we shall obtain that 



tf{X) =r]{XI-A,)-^BJ 



X+Wi 



6 



77 diagi 

<nS{Wn)* in 

E.ATfe-^.(^i-^(«^.)*^i'5(A)) 
Hiiir}iKs{X,Wi) 



(TiS 
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and consequently, 

{tf, tf) = ^ir]iKs{X, Wi),^. CjTiiKsiX, Wj)) = 

= (/,/)• 

Thus t is an isometry and Theorem 13.21 allows to conclude. 



□ 



The following result shows that the assumption X > in the statement 
of Corollary 13.31 can always be achieved for ri = 1. 

Lemma 3.4. Given S G TZS which is not the function identically 
equal to Ip. Then there exist a pair {^,w) G C^^^ x C+ such that the 
corresponding X > 0. 

Proof: We proceed by contradiction. Assume that for each w G C+ 
and for each vector ^, it holds that 

or, equivalently, 

CKs{w,w)C = 0. 

Then, ^^^(A,^;) = 0, and for each / G HiS) 

^f{w)=<f{X,^Ks{X,w) >=0. 

The space 'H(>S') is thus trivial, and its kernel is zero: 

a, - S{wya,S{X) _ 
X + w 

from where we conclude that for each w, X 

S{X) = a^^S-*{w)a^. 
and consequently, S{X) = Ip. □ 

3.2. Analysis of the tangential Schur algorithm. For n = 1, the 

matrix function in (13. lip becomes 
(3.12) 
0(A) 



hp- B*X-\XI - A)-^BJ 



h, 



2p 



-ax 

ai J X(A + w* 



X(A + w* 
"x(A + ty*) 



X(A + w*) 



X(A + u;*) . 



en 912 

©21 ©22 
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where using the Lyapunov equation (13.101) 

AX + XA* = X{-w* -w) = 

= -BJB* = -[-v e ] 
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' -ai " 




—T]* 


ai 







we find that 
(3.13) 



X 



w + w* 



To the best of our knowledge, Theorem 13.51 below is new and uses a 
realization theorem of M. Livsic |BLi] . It will be of much use in the 
following sections. 

Theorem 3.5. Let Sq E S with minimal realization 

(3.14) So{X) = Ip- B*X,\XI - Ao)''BoCT,, 

and let 9 be given by (I3.12p . Then S = Tq{So) G 7Z.S and a minimal 
realization of S is given by 

S{X) = Ip- B*sX-s\XI - As)-'Bsa,, 



where 
(3.15) 

(3.16) 



Bs 



Xc 



Bo 

v-^ _ 

Xo 
X 



Xo 




(3.17) 



A, 



Ao 

-r]aiB*X^^ 
Proof: From the definition 



-w 





^aiC - WiV* 
w + w* 

X 



S{X) =Te{So{X)) = 

= (Gn + So{X)Q2i)-\ei2 + So{X)Q22) 



rj r]ai 



5o(A): 



X 



X{X + w*) 'X(A + «;*)' 

X(A + ^*) ^"^"^^^^^^ X(A + ^*) 
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Let US denote here So{\)^* = a*; then the preceding expression be- 
comes 

r]*riai a*rj(Ti 



S{\) = (/,+ 



X(A + w*) X{\ + w*)' 



[t]* - a*)ria 



1^-1 



X(A + w*)' 



Using (/ + 6*a)-i = I-b*a 
S{X) = (/,- 



X(A + i(;*) 
1 



(^o(A) 



X(A + w* 



1 + a6* 
(?7* - a*)?7cri 1 



we get 



X(A + «;* 



770-1(77* - a* 



-)(^o(A) + 



(77* - a*)^eri. 
X(A + «;*) ■ 



5'o(A) 



X(A + m;*) 



X(A + w*) + 77^1(77* - 5'o(A)e) ' 

Suppose further that there is a reahzation of the form (I3.14p for the 
function 5*0 (A). Inserting it here we shall obtain 

c(\\ ^ c (\\ ~ 'S'o(A)r)[Wi^o(A) - g^i] _ 

^ ' °^ ^ X(A + 7/;*) + 77^1(77* -S'o(A)e) 

= Ip-B*Xo\XI-Ao)-'Boai- 

(77*-^o(A)r)[wi^o(A)-e^i] 

X(A + w*) + 77^1(77* - So{X)^*) ■ 

Let us denote 

M =X{X + w*) + r]ai{7]*~So{X)e) = 

= X(A + w*) + 77^1(77* - C) + B*oX,\XI - Ao)-'Boa,C- 

The preceding formula for S'(A) becomes 

(77*-5o(A)r)[wi^o(A)-e^i] 



SiX) = Ip-B*X,\XI-Ao)-'Boa, 



Xn^ 



^0 





X(A + w* 

X 



1 



-5a~^ 1 



77^1(77* -5o(A)e) 



Bo 
77-^ 



where 

a = XI — Aq 

[5 =-Boa,e 
6 =r]a^B*X^\ 

By definition of M we have 

M = X(A + w*) + 77(71(77* - r) - 5a-^/3, 
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and hence we obtain from the formula 
-I -1 



a P 
5 D 



1 



D - 6a-^/3 
that the last expression for S{X) is 



-6a~^ 1 



S{X) 



1 



X 



X,' 

M L ■' ^ J [ 1 

-5a~^ 1 



X 



X, 



-1 











X 



a 



/5 





-1 











A/ -A 



X 






X 



0"! 



X 



X 



-1 



B, 



ai. 



Thus we have obtained the realization fl3.15p - fl3.17p . Furthermore, the 
Lyapunov equation fl3.10p holds since 



AsXs + XsA*s 

An U 

+ 



X 



Xo 
X 



+ 



Xo 
X 



X 



X ^ _ 



X 

AqXo + XoA* -Bo(Ti{r]* -C) 



-Bsa.B^ 



s- 



The last equality follows easily from the Lyapunov equation for the 
given realization of and the formula fl3.13p . □ 
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4. The tangential Schur algorithm in the class SX 

Our strategy to the tangential Schur algorithm in the class SX relies 
on the following theorem. This theorem shows in particular that one 
cannot use the naive approach of applying the classical tangential Schur 
algorithm for each (that is, looking at as a mere parameter). 

Theorem 4.1. Let us fix the parameters Ox^a^, and 7, and the interval 
I. Then for every ^2 ^ I there is a one-to-one correspondence between 
pairs (7*, 5") such that S G SX and 7* continuous in a neighborhood of 
t2, and functions Y{X) G S{t2). 

Proof: Let (p be the map which to a pair (7^,, S) associates the func- 
tion S{\,t2) G S{t2)- The converse map '■ ^(•^) ~^ '^i^ih) was 
introduced in [M], [MVTl chapter 7] in a more general setting, and is 
defined as follows. Suppose that we have realized the transfer function 
Y{X) in the form 

Y{X) = - B*X^\XI - A,)-'Boa,{tl). 

Then we construct -8(^2) from the differential equation with the spec- 
tral matrix parameter Ai: 

(4.1) ■^[B{t2)ax{t2)]+A,B{t2)a2{t2)+B{t2h{t2) = 0, S(t°) = Bq. 

0,12 

In fact, the function B{t2) is given by the formula 

(4.2) B{t2) = - Ai)~^B^a^^^-\\ ta, ^2)^^^. 

Next we construct X(t2) on the maximal interval /, where it is invertible 
fl2.14p via the formula 

(4.3) 4-X(t2) = B{t2)a2{t2)B{t2)\ X(t°) = Xq. 
dt2 

Finally, we define 

7*(^2) = 7(^2) + a2{t2)B{t2yX-\t2)B{t2)a,{t2)- 
-cr,{t2)B{t2yX-\t2)B{t2)a2{t2). 

Then easy computations show that the function 

(4.4) S{X,t2) = Ip- B{t2rX-\t2){XI - Ai)-'B{t2)a,{t2) 

is in the class CX corresponding to the parameters o"i(t2),o"2(^2),7(^2), 
and 7,(^2) for ^2 e I. ????????????????? Moreover, the Lya- 
punov equation (13.81) holds for every t2 G I, and thus the real- 
ization (14.41) is minimal for every t2 G I. 



ON THE CLASS SI 



21 



Note that the composition (p o ip = id^ since starting from a function 
Y G Sit^), constructing Y{\,t2) and taking its value at we shall 
obtain again Y from the initial conditions of the differential equations 
flO) and (gJl) defining B{t2),X{t2). 

In order to show that tpocf) = id, we start from a function S'(A, 12) € «SX 
and take its value S{X,t2) G «S(t2)- Using the construction above, we 
shall obtain a function Y{\,t2). Note that the two functions S{\,t2) 
and Y{\,t2) have the same value at and maps solutions of the same 
input LDE to (possibly different) output LDEs, i.e. : 

s{\M) = <^.{\t2A)s{\ti)^~\\t2A). 
y(A, t2) = $UA, t2, 4)s{x, tl)<^-\x, t2, ti). 

Then the function ^"^(A, t2)^(A, ^2) is equal to Lp at infinity and is 
entire. By Liouville's theorem it is a constant function and is equal to 
Ip. Thus 

<l>,(A,t2,t^) = <(A,t2,t^), 

from where we obtain that 

<i>:\x,t2,f2)Ki>^,t2,4) = ip. 

Differentiating both sides of this last equation we get to 

=^[$;HA,t2,t^)<fUA,t2,t^)] = 

= ^'\X,t2,4)a^\t2)i-^it2)+iit2))Ki>^,t2,4). 

Since the matrices $*(A, t2, ^2)) ^'*{^^ ^2)? (^lih) are invertible we ob- 
tain that 7(t2) = l'{t2). □ 

Remark: Last theorems claims that the correspondence is between 
"initial" values S'(A,t2) and pairs ('S'(A, ^2), 7*(^2))- Notice that it is 
possible to obtain functions with the same 7*(t2) with different initial 
conditions: 

Proposition 4.2. Suppose that there exists a function Y G S{t2), 
which commutes with $(A,t2;^2) ^'^^ suppose that a function S G SI 
corresponds to certain vessel parameters o"!, a"2, 7, 7*. Then the func- 
tion SY belongs to the class SX and corresponds to the same vessel 
parameters o"i , a"2, 7, 7* . 

Proof: Using formula 12.211 we obtain that 

S{X, t2) = $.(A, t2, tl)S{X, tl)^-\X, t2, tl). 
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Consequently, 

SiX,t2)YiX) = $,(A,t2,t^)5(A,t°)y(A)$-i(A,t2,t^) 

intertwines solutions of the input fl2.16p and the output (12.171) ODEs 
with the spectral parameter A, and is identity at infinity, because S 
and Y and their product are such. Thus by the definition the function 
SY G ST and corresponds to the same spectral parameters as S. □ 
For example (to be studied in subsection 15. 5p . taking Sturm-Liouville 
(SL) vessel parameters 





' 


1 " 




'10" 




"00" 




1 





,0-2 = 





,7 = 


i 



we shall obtain that 



V 



-k{t2-tl) 







gfc{t2-iES) 



V- 



where 



V 



A Va 
1 1 



i 
X 



i.e. of the 



Taking "K, which commutes with cr^ ^(c"2A + 7) 
form 

2C(A) 

(4.5) Y{X) = h- ''^^> 

. c(A) a(A) 

for functions a(A), c(A), which are zero at infinity, we shall obtain that 
for any S G <SX, the function SY G SI and corresponds to the same 
vessel parameters. But in this case (which can be also generalized) it 
turns to be a necessary condition as the following theorem states: 

Lemma 4.3. Given SL vessel parameters cti, cr2, 7, 7*(2;), there exists 
a unique initial value S{0, A) up to a scalar t2-independent symmetric 
and identity at infinity function. 

Proof: For this it is enough to prove that if a function Y{X,t2) G 
SI intertwines solutions of the input LDE (I2.16P with itself, then 
Y{X,t2) = Y{X,t2), i.e is a constant function, which moreover com- 
mutes with $/(A, t2, ^2)- Indeed, if we are given two functions 5'i(A, ^2), 
5'2(A,t2), then the function Si^{X,t2)S2{X,t2) will intertwine solutions 
of the input LDE with itself and as a result is constant Y{X) G «S, which 
is moreover commutes with ^i{X, ^2)- Following the paragraph preced- 
ing the theorem, ^(A) = / and this means that ^i(A,t?) = ^2(A,t^). 
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Let us show first that if 5'(A, ^2) intertwines solutions of the input LDE 
fl2.16p with itself and is identity at A = 00, then it is a constant matrix, 
which commutes with $(A,t25^2)- Performing simple calculations, we 
can find that 



c^i ^(cr2A+7) 
where 



i 
A 



V 











Consequently, for the expression 5" (A, ^2) = *^'(A, ^2, ^2)'5'o(-^)'^^^^(''^) ^2, i'^^) 
to be identity for A = 00, it is necessary to "cancel" the essential sin- 
gularity arising from two entire functions $(A,t2,^2) ^^'^ '^'^^(A, ^2, ^2) 
(or more precisely, let them cancel each other). 
Using the formula for $(A,t2,^2) 

k{t2-t9.) 



^X,h,tl)So{X,tl)^-\X,t2,4) 



V 









e*=(*2-t§) 



v-'s{x,tt)v-' 









-k{t2-4) 



and considering coefficients of the exponents, we may conclude that in 
order to cancel the exponents it is necessary to demand 



v-'s{x,tl)v 



6(A) 
d{X) 



for some analytic in -\/A at infinity functions b{X),d{X) which are 1 
there. From here it follows that 



S{X,tl) 



(6(A) + d(A)) -(ci(A)-6(A)) 



A 



d{X) - 6(A) 



iKX) + d{X)) 



denoting here 1 — a(A) ~ ~ y 'j^iKX) + d{X)) and — c(A) = d{X) — 6(A), 

we shall obtain that 6*0 (A) is of the form (14.51) . i.e. commutes with the 
fundamental matrix $(A,t2,^2)- 

Finally, let us consider a function of this form and will require the 
symmetry condition for it: 



S*iX,4)a^SiX,4) = a,. 
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Plugging here the expression (14. 5p for S{X,t' 



a{\) 



ic{\) 
X 



c(A) a(A) 



-A, the 



we shall obtain that on the one imaginary axis, where A = 
following system of equations must hold 

(1 - a*(-A))c(A) + c*(-A)(l - a(A)) = 0, 

(1 - a*(-A))(l - a(A)) + ^c*(-A)c(A) = 1, 

f 

(1 - a*(-A))(l - a(A)) + ^c*(-A)c(A) = 1, 
zc*(-A)(l-a(A)) ^ zc(A)(l-a*(-A)) ^ ^ 
A A 
Subtracting the second and the third equation we obtain that (remem- 
ber that A = —A) 



-A)c(A)(i4) 



;-A)c(A)- 



0. 



from where it follows that on the imaginary axis at the points of an- 
alyticity of c(A) it holds that c(A) = 0. Since it is analytic there, it 
will be globally zero. Notice that the last system of equations for the 
trivial c(A) gives 

(l-a*(-A))(l-a(A)) = l, 

Which means that 1 — a(A) is of norm one on the imaginary axis. □ 
We now focus on the rational case. We first note that starting from dif- 
ferent realizations at we shall obtain the same 7*(t2)- More precisely 
we have the following theorem: 

Theorem 4.4. Suppose that there are two minimal realizations of the 
function S{X,t2) E IZS: 

S{X,tl) = Ip- B;X-\XI - A,)-^Beai{tl), £=1,2, 

with associated similarity matrix V. Then the functions 

Si{X, t2) = Ip- Be{t2rXi\t2){XI - Ae)-'Be{t2)a,{t2), £=1,2, 



obtained from these realizations via the construction in Theorem \4.1 
are the same, and this holds if and only if: 

(4.6) B2it2) = VBi{t2), 

where t2 varies in a neighborhood of Moreover, in this case the 
following formula holds 

X2{t2) = VXi{t2)V*. 
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Proof: Equality of the two realizations means that there exists an 
invertible matrix V such that 

A2 = VAiV-\ B2ai{tl) = VBiaiitl), BIX^W = B;X^\ 

from which follows (see \AG\ Lemma 2.1 p. 184] for instance) that 

X2 = VXiV*. 

By the construction described in Theorem 14. 11 the function Bi{t2) will 
satisfy (jH]) 

-^[Si(t2)ai(t2)] + A,B,{t2)a2{t2) + B.it^Hh) = 0, B,{4) = B,, 
at2 

and the function -82(^2) will satisfy the same equation with A2 instead 
of A: 

^[B2it2)a,it2)] + ^252(^2)^2(^2) + B2{t2hit2) = 0, 52(t°) = B2. 

at2 

Using the equalities A2 = VAiV~^, B2 = VBi we obtain that the 
function ^"^^2(^2) satisfies 

^[■^^-'^2(^2)^1(^2)] + AiV-'B2{t2)a2{t2) + V-'B2{t2h{t2) = 0, 
V-'B2{t'i) = 

which is the same differential equation as for Bi(t2). Thus ^2(^2) = 
VBi{t2)- Similarly, considering the differential equations 

4rMt2) = B,{t2)aiB,{t2)\ X,(tO) = X„ z = 1, 2, 
dt2 

we obtain that 

X2(t2) = VXi{t2)V*. 

Consequently, 

S*(t2)X2 '(t2)5(t2) = Bl{t2)V*V-^*X^\t2)V-^VBi{t2) 
= Bl{t2)X^\t2)B^{t2), 

from where we conclude that the same function 7^. is associated to 
Si{\,t2) and S'2(A,t2)- Since these functions coincide for t2 = ^2 
map the same input ODE, we obtain that they are equal in a neigh- 
borhood of ^2- n 

The following notion has been introduced and studied in [M2j. 
Definition 4.5. Let S G IZSX with a realization fl4.4p . The function 

T{t2) = detX(t2). 
is called the t -function associated to S . 
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It follows from Theorem 14.41 that the r function is well defined up to a 
multiplicative strictly positive constant. Indeed using the notation of 
the theorem, 

detX2(t2) = det[VXi{t2)V*] = detXi{t2) det{VV*). 

We now introduce the counterpart of the tangential Schur algorithm 
in the class IZSI. Let 5* G IZSI and fix ^2 ^ I- The function 
S{X, t^) £ TZS. Consider now a space with X > and corresponding 
function B as in Corollary [331 This is always possible in view of Lemma 
13.41 It follows from Theorem 13.21 that there exists 5*0 G IZS such that 

(4.7) S{X,4)=Tew{SoW). 

Applying Theorem 14.11 to 5*0, we obtain a uniquely defined function 
So{X,t2) G IZSI, such that at the relation (14. 7p holds. 

Definition 4.6. The map T^ q 

is the time-varying counterpart of the linear fractional transformation 
(14.71) . We will call it a generalized linear fractional transformation. 

If S{X,t2) G SI corresponds to the vessel parameters cti, a2, 7, 7*, 
then So{X,t2) G SX corresponds to the vessel parameters ai, 02, 7, 
7o,* for a uniquely defined function 70,* (^2)- Moreover, for t2 = we 
have the usual linear fractional transformation (14.71) . 

As a consequence of Theorem 14.11 the following lemma holds. 

Lemma 4.7. For a given J -inner function and a given point t2, the 
map Tq ^o is one-to-one from IZS into IZSI. 

Proof: Notice that for a given ^2 the map Tq is injective. Further- 
more, using Theorem 14.11 every S G IZSI is uniquely defined by the 
function S{X,t2) G 7ZS. The result follows. □ 

Suppose now that we start from So{X) = Ip and apply n linear frac- 
tional transformations for a fixed t2, using the data < Wi,^i,r]i > 
{i = l,...n) to construct the corresponding J-inner functions. We 
obtain a function 

Sn{X) = Te(/) = / - B:X-\XI - A^y'B^a, G TZS. 
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Using iteratively formulas fl3.15p . f l3.16p . f l3.17p we obtain that 

Bn = : 
X„ = diag[Xi, . . . ,X„], 
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and 

(4.8) 



m^ii-nl -il) 



r?20-i(r)* - c*) 



iV2 - 



-W„ = 



where Xj is defined by (13.13P : 



X,- 



Wi + w* 



1, . . . Ti. 



Assume now that, starting from the identity matrix we apply this pro- 
cedure for two different sets of data (with the same n) at two different 
points tl and t^- The following theorem answers the question as when 
we obtain the same function, that is, when do we have: 



Theorem 4.8. Suppose that there are given two sets of n triples 
< w},^l,r]l > and < Wi,^f,rif >, with corresponding 6^,£ = 1,2. 
Then necessary and sufficient conditions for equality of the two func- 
tions 



1,2 



are: 



(1) The corresponding matrices A\ and A\ defined by (14.81) are sim- 
ilar, i.e. there exists an invertible matrix V such that A]^ = 
VAlV-\ 

(2) §{\I - A\)-'Bia-,'^-\\tlA)dX = VBl. 



Proof: From Theorem 14. H a necessary and sufficient condition for the 
functions to be equal is that 

s,{\tl) = S2{\A)- 

From 14.41 this holds if and only if 

Al = VAiV-\ Blitl) = VBl Xlitl) = VXlV* 
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for a uniquely defined invertible matrix V. The result follows using 
formula (jtg). □ 
A more general construction in this setting is obtained if one supposes 
that at each step different values of t2 are chosen. In this case the 
construction of the function S'„(A,t2) is more comphcated, and can be 
computed recursively. The formulas are very involved in this case and 
we can see no real advantage to develop them at this point. 



5. Markov moments and partial realization problem in 

the class sx 

Let S E S, and consider the Laurent expansion at infinity 

oo ^ 

S{X) = Ip- B*X-\XI - A)-^Bai = Ip-J2 -^B*X-^A'Bai 

in terms of a given minimal realization. The matrices Hi = B*lL~^A'^Bai 
are called the Markov moments of S. The Partial Realization Prob- 
lem (or moment problem at infinity) in S is defined as follows: Given 
the first n + 1 Markov moments Hq, . . . , Hn, find all functions (if any) 
5" e «S with these first n+1 moments. See |GKL] for a general study of 
the partial realization problem. Similarly, one can define the Markov 
moments for an element S G SX. We now give necessary conditions 
which the moments of a function 5* G SX have to satisfy. 
It is important to notice the following: fixing t2 = ^2 ^"^^ solving the 
corresponding classical moment problem will not lead to a solution of 
the problem in the class SX because we cannot obtain the function 7 
(and hence 7*) from this solution. 

5.L Fundamental properties of S{X,t2). Before we consider mo- 
ments of the function S{X, ^2) we present some of its fundamental prop- 
erties which shed more light on the deriving of the moment equations. 

Theorem 5.1. For fixed ai,o"2,7, a necessary condition on 7^, so that 
there exists a vessel with the vessel parameters 0\,02, 7, 7* is 

det ^^(A, t2, ^2) = det <I>(A, ^2, ^2) 

Proof: Let S G SX be a function corresponding to the parameters 
cr2 5 7, 7*- Then using a realization 04.41) 

S{X, t2) = Ip- B{t2yX-\t2){XI - A,r'B{t2)a,{t2) 
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and Lyapunov equation (13. 8p we shall obtain that 

det S{\M) = det (Ip - S(t2)*X-i(t2)(AJ - Ai)-^B{t2)ai{t2)) 
= det (/ - B{t2)aiit2)B{t2)*X~\t2){XI - A^)-^) 
= det (/ + {AiX{t2) + X{t2)Al)X-\t2){\I - A,)-') 
= det (/ + Ai{XI - Ai)-^ + X{t2)AlX-\t2){XI - A^y^) 
= det ((A/ + X{t2)AlX-\t2))iXI - 
= det(A/ + Al) det(A/ - A^)-^ 
= det S{X,tl). 

Consequently, taking determinant of the formula fl2.21l) 

six,t2Mx,t2,4) = <^.ix,t2,4)six,4) 

we shall obtain that det $*(A, t2, ^2) = det $(A, ^2, ^2) for all points A, 
where det S'(A, ^2) exists and is different from zero. Since it happens for 
all points outside the spectrum of Ai and the functions det $*(A, t2, t^); 
det $(A, t2, ^2) are entire they are equal for all A. □ 
Next theorem put some light on the connection between equations, 
which determine the transfer function S{X,t2) and its moments. This 
theorem is similar to the property of a solution of a Riccati equation 
\Ze\ theorem 2.1] 

Theorem 5.2. Suppose that S{X, ^2) is a continuous function of t2 
for each X, meromorphic in X for each t2 and satisfies S{oo,t2) = I. 
Suppose also that S{X, ^2) is an intertwining function of LDEs (12.161) 
and (I2.17P . Then if the symmetry condition (12.231) 

S{XM) = ^i\t2)s-^*{-x,t2)(yi{t2) 

holds for t^, then it holds for all t2- 

Proof: Since S'(A,t2) intertwines solutions of (12.161) and (12.171) . then 
it satisfies the differential equation (12.221) 

J-^(A,t2) =CTfl(t2)(a2(t2)A + 7*(t2))5(A,t2) - 
012 

-5(A, t2)ay{t2){a2{t2)X + 7(t2)). 

Consequently, using properties of 7*, 7 appearing in definition 11.11 we 
obtain that the function crjf ^(t2)5'^^*(— A, t2)ci(^2) satisfies the same 
differential equation. If these two functions are equal at ^25 from the 
uniqueness of solution for a differential equation with continuous coef- 
ficients, they are also equal for all ^2- D 
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5.2. Restrictions on Markov moments for functions in SX. We 

study the Markov moments of a function S G «SX, which maps solutions 
of the input ODE (12.161) to the output ODE (12.171) using a minimal 
realization (14. 4p of 5". 

If Ax is a constant matrix, at some stage the elements /, Ai, . . . , ^4" 
will be linearly dependent and we obtain: 

Lemma 5.3. Given S G IZSI, with Markov moments Hi{t2), i = 
0, 1, . . .. Then exists N and constants fij, j = 1, . . . , N such that 

N+l 

(5.1) ^'lHn-,{t2) = 0, n > + 1. 

Next we present formulas, which hold for the moments of a function in 
SX. They can be easily derived from theorem 15.21 We notice that the 
first moment i?o(^2) satisfies the linkage condition (II. 3p : 

7*(^2) - 7(^2) = (r2{t2)Ho{t2) - cri{t2)Ho{t2)(rx\t2)(r2{t2). 
Let us denote 

Ho{t2) = C{t2)B{t2)a,{t2) = {B{t2)yX~\t2)B{t2)ai{t2). 
The functions B{t2),C{t2) satisfy the following differential equations 

(5.2) A[5(t2)ai(t2)] + AxB{t2)a2it2) + 5(^2)7(^2) = 

Clt2 

(5.3) ax{t2)^C{t2) - a2{t2)C{t2)A, - 7*(t2)C(t2) = 0, 

at2 

see jMl IM VI] . Thus, differentiating i/o(^2), we obtain 

^//o(t2) = ^^[C{t2)B{t2)a,{t2)] 

= (Tl\t2)a2{t2)C{t2)AxB{t2)ai{t2) - C{t2)AxB{t2)a2{t2) + 

+(ri\t2h,{t2)Ho{t2) - Ho{t2)cr^\t2h{t2). 

In other words the second moment Hi(t2) = C{t2)AiB{t2)ai{t2) satis- 
fies the following differential equation 

a^\t2)a2{t2)Hi{t2) - Hx{t2)ax\t2)a2{t2) = 

= i;Ho{t2) - ax'l.Ho{t2) + Ho{t2)ax\t2)i{t2). 

In the same manner the moments Hi{t2) and Hi^i{t2) are connected 
by the following differential equation 

(Ji\t2)a2{t2)H,+x - H,+i{t2hi\t2)a2{t2) = 
^^•^^ = £H,{t2) - a^\t2h,{t2)Hi{t2) + H,{t2)ax\t2)i{t2). 
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Notice also that 



Hia^^H*ai = B*X-^A\BaiB*X-^Bai = 

= B*X-^A{[-AiX - XAl]X-^Bai = 

= -Hi+i - B*X-^A{XAlX-^Bai = 

= -Hi+i - B*X-^A\-^[-BaiB* - XAl]AlX-^B(Ji = 

= -Hi+i + B*X-^A\-^BaiB*AlX-^Bai + 

+B*X-^A\-^XAlAlX-^Bai = 

= -Hi+i + Hi^ia^^H*ai + B* X'^ A^^ XAf X'^ Bai 
Consequently, the following formula holds: 



(5.5) i/,+iar^ + i-iya^'Hl, = Y,{-iy^'H,_,a^^H*^ . 

j=0 

Remark: One can obtain the same equations (15.41) and (15.51) by taking 
the expansion of ^(A, t2) into Taylor series around A = oo and equating 



the coefficients of each — . More precisely, the algebraic condition (jSl 

is a consequence of the symmetry condition (12.231) and the condition 
(15. 5p is a result of the differential equation (I2.22p . 
Finally, we show how the third condition in Proposition 12.51 is reflected 
in the moments Hi{t2). This condition means that the function S'(A, ^2) 
is (Ji{t2) contractive, and, for example, the first moment i/o(^2) satisfies 
Hoit2)(J^^(t2) > 0. Using the minimality property, we have 

span A'lBS = H. 

Thus, in order to have X > eJ, e > (so that X^^ exists) it is enough 
to demand 

n n 

(X^^ ^ A^Bek, ^ A^Bek) > e, for any choice of G £. 

k=0 k=0 

Thus we obtain that the following Pick matrix arises 

B* 



B*Al 



S*(A*)" 



X"^ [ B AiB . . . A^B ] > el 



which is equal to 



n 



[B*{AiyX-^A{B] > el, 



n 



0,1,2,... 
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Notice that the elements of the matrix are indexed from zero to n for 
convenience. Using Lyapunov equation (I3.8p iteratively, we shall obtain 
that 

m 

= E (-l)'+Vr'//:,_i_,aiif„+fcar\ m > O, 
B*X-^AIB = Hn. 

So the Pick matrix is of the following form 
(5.6) 



k=0 

PO,j =Hj. 



Next theorem shows that moments of a transfer function in S satisfy 
recursive equations, involving algebraic and differential equations: The 
positivity of the matrix P, in that case does not necessarily hold. 

Lemma 5.4. Denote the real and the imaginary parts of H2n+iO'i^ = 
R2n+i + 'iM2n+i- Then if the moments Hq, . . . , H2n satisfy the differ- 
ential equations (15. 4p and (15. 5p . then the real part R2n+i satisfies the 
imaginary part of differential equation (15.51) : 

di \ ~\ \ \ \ ]_ ]_ 

d 
dt2 



Proof: Multiplying equation (15. 4 p by ^ on the right and using the 
notation H2n+iO'i^ = i?2n+i + ^^2n+i we shall obtain that 

erf V2(i?2n+1 + iM2n+l) ~ (i?2n+l + «M2„+i)(T2Crf ^ = 

= af V2i?2n+1 - -R2n+lCr20-f ^ + i[o 2^^^ M2n+\ - M2„+l(T2af ^ ) = 

= —\H2y^(y^^ - (j];^-i^H2n<yi^ + i/2nO-f So"r^ 

CLt2 
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Since V2/22n+i - R2n+io-20-i^ and cr2crf ^M2„+i - M2n+io-2crf ^ are 
skew adjoint, last equation is equivalent to equality of real and imagi- 
nary parts of both sides: 
(5.8) 

2z[(Tf V2M2„+i - M2„+ia2fTf 1] = —[H2n]a^^ + (J^^ — [H*J + 



and flO) 



2[ai V2i?2n+1 - -R2n+lO-2Cri ^] = Tr[-^2n]o"i ^ - CT^ S*-f^2nC^l ^ + 

(112 



2n- 



Let us show that i?2n+i, defined in equation flS.Sp also satisfies (15. 7p . 
Indeed, 



2[cri V2i?2n+1 — -R2rt+lO"20"i ^] = 

= —[H2n]o-i^ - Cr^^J*H2nO-i^ + H2n(yi^ ICf^^ - 
0,12 



Performing rearrangement and using flS.Sp for i = 2n — 1, we shall 
obtain 



^i'^2n:T=oi-^y^'H2n-,a^'H*] - [ET=oi-^y^'H2n-,a^'H*]a2a^' = 

= —[H2nai^ - (^l^H^J - a^^'J^H2n<Ji^ + H2nCri^l<7i^ " i^2nTr'^r^" 
(112 0,12 

= ^^[j:T=o\-^y^'H2n-l-J^l'H*] - a^'j,H2na^' + i/2n^rS^r'- 
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Here we can use formula (15.41) in order to differentiate the expression 

H2n-i-jO-i^H*. From ( 15^ it follows that 

(5.9) 

= [cTi Vai/fe+i - Hk+ia^ Vs + ^l^Hu - H^a^ ^'j]a^ ^H^+ 

= cr2-fffe+icrf + HkCTi H^^^a2<Ji ~ 
—Hk^ia^^a2(Ji^H*^ — iffecrf V20-f + 

Now we can insert the formula (15.91) into the left hand side of the 
previous expression and we shall obtain 



= ET=o\-^y^'[^i'^2H2n-,cr^'H; + H2n-i-jCr^'H*^'a2a^'- 
+a^'j,H2n^i-j(J^'H* + H2n-i-j(Ti'H*i:a{']- 



It is easy to recognize telescopic sums in this expression and we obtain 
that finally it equals to 

= ET=o\-'^y^'l^i'^^H2n-j(r^'H* + H2n-i-ja^'H;^,a2a^'] + 
+H2nCX^^cr2a^^H* - Hoa^'a2a^' H*„ + 

= ET=oi~^y^'[^i'^2H2n-,a^'H; - H2n-jaY'H;a2a^S + 

+ Epo\-^y^'Wi'^*H2n-i-ja^'H; + H2n^i-,a^'H*i:a^']- 

In other words, we obtain that the compatibility condition (15. 7p be- 
comes 

Y.T=o\-^y^'[^]'l*H2n-i-,a^'H* + H2n-i-,o^'H*i:a^']- 
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or after multiplying by ai from both sides 

2n-l 

□ 

Theorem 5.5. Suppose that we are given moments Hi{t2), defined in 
a neighborhood of the point ^2 ^ I. Then there exists Uq < such that 
each element Hi^i is uniquely determined from HQ,...,Hi using the 
algebraic formulas fl5.5p . andriQ LDEs with arbitrary initial conditions, 
obtained from (15.41) . Moreover, equations must be satisfied by the 
elements o/7*(t2)- 

Remarks: 1. This theorem is of local nature. The number "no" ap- 
pearing in the theorem may vary with the point t\, but is unchanged 
in a small neighborhood of t^ by continuity. 

2. The proof of the theorem allows to produce an algorithm to deter- 
mine explicitly, up to no initial conditions, the Markov moments. The 
arguments in the proof of the theorem are illustrated on an example 
in the following subsection. This example exhibits all the difficulties 
present in the general case. 

Proof of l5.5t We have seen in lemma [5^ that the real part of H2n+i = 
R + iM, defined by equation (15.41) satisfies (15.71) . Let as concentrate 
on the equation (15.81) in order to solve for M. For example, if 

spec <yi^a2 n spec(— cr20"f ^) = 

then one can uniquely solve for M. But in general, there will be some 
undefined entries of M in this manner. In the later case, we can rewrite 
these equations as linear equations in the entries of M. But before we 
do that we also rewrite the equation (15. 8p as following 

(5.10) 2z[arV2Af2n+l - M2n+l(J2(Tl'] = 

= -TrR2n - (^l^l*H2nCri^ " H2nO-i^ J* CT^^ - 
0X2 

And further, we are able to rewrite this equation as a system of 
linear equations in variables M^', k,l = 1, . . . ,p'^: 



(5.11) 



az2 



'2nJ 



[ E c^^MU = i-RZ + YXfiZR'in] + Y.WiMt] 
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Here a^i, 6]^i are functions of t2 derived from the vessel parameters. 
Using basic algebra manipulations, it is enough to find a maximally 
independent subset of equations in the left side of (15. lip , and to obtain 
a system of — hq independent equations 



(5.12) 



f ( ^ T^kl TDkl i\/rkl\ 
J2h '^fci (^2)^/2^^+1 = /2( — -^2!) M2n) 



2n+l 



/p2-no(^^ -^2rn -^2ni ^2n) 



for linear functions and a system of linear dependent ones 



d 



— 9li-^R2n, R2n^ ^"2^ 
= g2{-J^R2n, R2L ^2n) 







,_d_ 

'dt2 



9no [ — R2n, R2n, M^'' 



2n) 



for linear functions Qj. 

On the other hand, using the same considerations for H2n+2, "we shall 
obtain hq equations of the same type, which are additional restrictions 

on M2n+1 



(5.13) 



d 

= 9j{-rrM2n+l,M2n+l, R2n+l)J 
0,12 



.no. 



Notice that we obtain algebraic (I5.12p and differential (I5.13P equations 
which are independent (since all algebraic and differentials are, and 
they are of different nature). So, if all elements of M2„+i appear at 
these equations we shall obtain that one can uniquely solve them up 
to no initial conditions: 



(5.14) 



d 



kl\ 
2n) 



d 

= 5'i(Tr^2n+l, ^2n+l, -R2n+l), 
dt2 



,P - no, 



1, . . . ,no 



(See example of NLS at section [5l6!) . If this is not the case, some of 
the elements, say po < no will not appear in the equations (I5.14p . As 
a result, using the algebraic expressions (I5.12p and plugging them into 
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differential ones, f l5.13p we shall obtain po < differential equation of 
second order on the elements i?2n and M2n'- 



d 



' dt2 dt2 



dU 



By induction the same type of equation will hold for if2n+i (for j 
l,...,Po): 



d 

dt2 



d^ ^ 

(5.15) = hi{-^M2n+l, —M2n+l,M2n+l, — i?2n+l , ^2n+l) • 



d 

dt2 



Finally, we shall obtain the following system of differential equations 



' E..«L(^2)M|^+i = M^^Rll,Rll,Mt) 

^ 2 = 1, . . . - no, 

(5 16) '\ *^ ^ fi'i(^^2n+l, A^2n+1, R2n+l) , 

1 = l,...,no-po 
d^ d d 

^ ~ ^j("7;:2^2n+l) -7;-^2n+l) ^2n+l) -7;--R2n+l) R2n+l) 1 
0X2 dt2 (lt2 

2 = 1, . . . ,P0 



□ 

Remark: This theorem means that the equations (15 ■4p and (15.51) are 
compatible, i.e. one can always solve them, once no equations are 
satisfied for 7*(t2)- But if we consider the differential equation (15. 4p 
only and ignore (15. 5p . then one obtain the following 



Corollary 5.6. Suppose that Moments of Hi satisfy the differential 
equations (15.40 . then there exzs^ -^"0(^2), -^1(^2), -^2(^2), -^3(1^2) depend- 
ing on vessel parameters only such that 

(5.17) ^^n+l = ^O^n + Kl^Hr. + ^2^^^n + ^3^^n 

Proof: All we have to do it to imitate the proof of theorem 15.51 ignor- 
ing the additional constrain (15. 5p on the real/imaginary parts of the 
moments. Denoting by H^l the entries of iJ^, we will obtain from (15. 4p 
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equations (like (15. lip ) 



n+l 



,12 TTkl 



n+l 



d 

d 

dU. 



d 

dti 



and analogously to deriving the formulas (15.141) . we will obtain 



d 



= gi{—Hn+i, Hn+l), 

0,12 



l,...,p2. 

1, . . . ,no 



^0, 



So, if these equations are independent, differentiating the first one, we 
shall obtain 

d d d"^ 

dt2 0X2 UX2 

Otherwise, there will be po equations analogous to (I5.15P of the form 

d^ d 
dt2 dt2 

Moreover, the independent entries of i^n+i will appear with one differ- 
entiation only, because, otherwise they would appear at the previous 

stage at the equation = gi{——Hn+i,Hn+i), and would be derived 

dt2 

from no differential equations. So, this last equation means that the 
remaining Pq elements are found from 



d 



d 



dt 



■dti 



L2 Ut-2 "''2 

SO that only the known entries appear at the left hand side. From here 
it immediately follows the corollary. □ 
We want to present next a necessary restriction on 7*(t2), derived from 
the existence of a finite dimensional vessel: 

Theorem 5.7. Let cri,cr2,7,7* be vessel parameters, and QJ a finite 
dimensional vessel corresponding to the m parameters. Then the en- 
tries of the function 7* satisfies a polynomial differential equation of 
finite order with coefficients in the differential ring IZ, generated by 
(the entries of) ai, ctj^"'^, a2, 7. 
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Proof: Suppose that the transfer function of the vessel 2J, defined in 

(En]) is 

oo 

S{X, h) = I- B*{t2)X-\t2){XI - A,)-'B{t2)a,{h) = I 

1=0 

and the hnkage condition is 

Cri\t2h4t2) = CT^\t2hit2) + [crY\t2)a2it2),Ho{t2)]. 

Notice that if we differentiate this formula and use the equation for the 
derivative of i?o(^2) from the equation (15 ■41) . we shall get 

^K^(t2)7*(t2)) = ^(ar^(t2)7(t2)) + [^^{a^\t2)a2{t2)),H,{t2)] + 
+ [a^\t2)a2{t2),^Ho{t2)] = 
= ^(^r'(t2)7(t2)) + foo{Ho{t2)) + foi{Hiit2)) 

for linear in moments functions /oo, /oi with coefficients depending on 
IZ and 7=K. Similarly, differentiating this expression and using formula 
d d 

(15.41) for "1— -f^o(^2) and ——Hi{t2), we shall obtain that the second 

uC2 dt2 

derivative is 

^(ar'(t2)7*(i2)) = 

= ^K'(t2)7(t2)) + flo{Ho) + fn{H,{t2)) + /l2(^2(t2)). 

for linear in the moments functions /io,/ii,/i2 with coefficients de- 
pending on IZ and 7*, ^—7*. Continuing this differentiation further at 

dt2 

each step i we shall obtain an equation of the form 

^{a^\t2h.{t2)) = ^{cr^\t2h{t2))+i2UH,): 
2 2 ^=0 

where fij is a linear function of Hj with coefficients, depending on TZ 
and first i — 1 derivatives of j^, (this can be immediately seen by the 
induction). But at some stage the moments start to repeat themselves 
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due to equation (15. ip . So, taking K derivatives we shall obtain equa- 
tions of the following form 
(5.18) 

f (yi\t2)l,{t2) = a^\t2h{t2) + foo{Ho{t2)) 

j=o 



AN jN N 



(it 2 (it 2 j=o 

fjK jK N 

;^K'(t2)7*(t2)) = ^K'(t2)7(t2)) + E //..(^.O 

j=o 

Suppose that that the dimension of the inner state is n (i.e. dimH = 
n), then we get that each of the matrices Hj has entries and there is 
the total number of n^(A^ + 1) entries for the moments Hq, . . . , H^. So, 
taking "enough" derivatives of 7=1,(^2) (i-e., taking K so that K dim{Sy > 
ii?{N -\- 1)) and eliminating all the entries of the moments, we shall ob- 
tain a finite number of polynomial differential equation for the entries 

of7,(f2). □ 
Remark: From this theorem is follows that 7=,, satisfies an equation of 
the form 

P(x,x',x",...,x(-^)) = 0, 

where P{xo, xi,X2, ■ ■ ■ , xk) is a non-commutative polynomial with co- 
efficients in 7?.. 

Using an outer space transformation it is possible to bring the matrix 
o"r^(^2)c"2(^2) to a simpler form. Suppose that there is an invertible 
matrix V{t2) : £ £ such that 

a^\t2)a2{t2) = V{t2)L{t2)V-\t2) 

for a "simple" matrix L{t2), where ^^(^2) is in a Jordan block form. 

Define new operators (omitting t2-dependence) 

(5.19) 

ai = V*aiV, ?2 = V*a2V, 7 = V*-fV - aiV, 7* = V*-f,V - a^V 

then simple calculations show that 51(^2), ^2(^2)5 7(^2), 7*(^2) are vessel 
parameters in the same interval I and the collection of operators and 
spaces 

^= (Ai,S(t2),X(t2); 5i(t2), 5^2(^2), 7(^2), 7*(^2);^,^), 
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is again a Vessel. By multiplying the original vessel conditions of 5J by 
V or on V* we check that the the vessel conditions are satisfied for QJ. 
But for this new vessel 5J we obtain that 

The first step to understand the equations arising in theorem (15. 5p we 
study the cosntant case. 

5.3. Moment equations for diagonal constant af ^cr2. At the spe- 
cial case o"f V2 = diag[si, S2, . . . , s^;] {E = dim£) for some constants 
Sj e C we can explicitly find the formulas for the moments using equa- 
tion <^;^. 

First we consider the linkage condition (11.31) . which may be considered 
as the " — 1" equation of (15.41) . Taking diagonal form for 0"f^(t2) 0^2(^2) 
and denoting the fc, j-th entry of i/o(^2) by Hq-', we shall obtain 

<y^\h)b.{h)-i{h)] = ^-'(^2)^2(^2), ^0(^2)] = [ H^'{t2){sk-s,) ] 

Cases for which s^. — Sj = give restrictions on the entries of 
c"r^(^2)[7*(^2) — 7(^2)] and other cases uniquely determine entries of 

Ho{t2y. 

(5.20) efcari(t2)[7*(t2) -7(^2)]e5 =0' - Sj = 0, 

ir, on Tjkn, \ efc^r'(*2)[7*(i2) - 7(^2)]e* 

(5.21) HQ-'{t2) = Sk-Sji-Q. 

Sk Sj 

But there more equations arise when we consider f l5.4p for i = 0: 

K"V2,//l(t2)] = [ H'^\h){Sk-S,) ] = 

ai2 

Here again if we focus on the cases Sk — sj = 0, we shall obtain that the 
following differential equations must hold (suppressing notation for t2 
dependence) 

^Hq^ - efc[o-]"S*^o + ^^ocrrS]e* =0, Sk - Sj = 0, 
^^^'(^2) = ekl^Ho - arS*^o + Hoa^^^U, Sk - Sj ^ 

Sk — Sj (112 

The first line gives a linear differential equation for the elements Hq-' 
in terms of vessel parameters and other known entries of Hq. 



-j-Hq^ = CkCF^ S*^oe* - ekH^a^ ^e*, 1 < A; < r, 

(112 



Finally we obtain the following 
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Lemma 5.8. Suppose that 

cp-Vs = diag[si,S2, . . . ,se], 

then the entries 0/7,7* satisfy the following equations 

(5.22) efca]"^(t2)[7*(^2) - l{h)]ei = 0, when Sk - Sj = 0, 

where stay for the standard elementary row vector with 1 at k-th 
place and everywhere. 

An interesting question arises, is whether any continuous vessel pa- 
rameters satisfying the conditions of this lemma gives rise to a transfer 
function. 



5.4. Moment equations for Jordan-block constant ^a2- Sup- 
pose now that ^o"2 is a Jordan block matrix with arbitrary eigenvalue 
s e C. 

s 1 



0"i V2 



Jordan(s, r) 



s 













s 1 
s 



then we can again consider the linkage condition (II. 3p : 

^r'(t2)[7*(i2) -7(^2)] = W^\t2)cX2{t2),Ho{t2)] = 

= [Jordan(s,r), 7/0(^2)] = [Jordan(0, r), i7o(t2)] 





K 


-Hi' 


TT23 

-"0 


- Ho 


• • Hl^ 


- H','-' 




HI' 


-Hi' 


TT33 


7T22 

- 


■ ■ Ht 




K 


K 


-Hf 


TTA3 

-"0 




■ ■ Ht 


- nr-' 



-1,1 






H, 







1,2 rrr-2,1 

— -Hn 



'0 

- H^o-''' 



-1,3 



^0' - H, 



Ho-'' 
r-1,2 



-l,r 



2,r-l 



-Ho 

- Hq 



where the value s on the main diagonal is canceled. It turns out that 
sums of elements of each diagonal under the main one is 0. This may 
be rewritten also in the following form 

tr[af i(7,(t2) - 7(^2)) Jordan(0, r)^] = 0, = 0, 1, . . . , r - 1. 

From the other equations we find expressions for the entries of -^0(^2) 
using its first row as parameters. In other words, from the linkage 
condition (11. 3p we find r restrictions on erf ^(7*(t2) — 7(^2)) and all the 
entries of i^o(^2), while the entries of the first row H^ 
unknown. 



11 
; 



Ir 



are 
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Considering next the equation (15 .4^ . we will similarly obtain that 



tr[(-^i7o - (7 1 ^1*Hq + H^a^ ^7) Jordan(0, r) 
at2 



0,fc = 0,l,...,r-l. 



These equations may be considered as additional equations on the first 
row of ifo(^2); which was unknown before. It is also possible that these 
differential equations will not determine the first row of Hq, in this 
case we shall obtain additional differential restrictions on the elements 
of (Tr^(7,(t2) - 7(^2)) (see (ICTll below). In this case, the elements of 
the first row of Hq will be determined from differential equations on 
the next level, involving Hi{t2) (see (15.321) below). 

5.5. Sturm Liouville vessel parameters. The following example 
was extensively studied in jM2] . It deals with the Sturm Liouville 
differential equation 



dtl 



with the spectral parameter A. The parameter ^(^2) is usually called 
the potential. For ^(^2) = this problem is easily solved by exponents 
and in this case we shall call this equation trivial. In |M2] one connects 
solutions of the more general problem with non trivial (7(^2) to the 
trivial one. 

Definition 5.9. the Sturm Liouville vessel parameters are given by 



7*(^2 



1 

1 



0-2 



1 




,7 



-27rn(t2) -/3(t2) 

/3{t2) I 





i 



for real valued continuous functions 7rii(t2), l3{t2)- 

Notice that these parameters correspond to the Jordan block form of 



0^1 cr2 





1 



The input compatibility differential equation (12.161) is equivalent to 



XUi{\,t2) - ^M2(A,t 



_d_ 

dt2 







where we denote U\{t2) 
one finds that M2(A,t2) 



-Ui(A,t2) + iu2{X,t2) = 

From the second equation 



Ui{X,t2) 
U2{X,t2) 



—i-£^ui{X,t2) and plugging it back to the 
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first equation, we shall obtain the trivial Sturm Liouville differential 
equation with the spectral parameter iX for Ui{X,t2): 



dt 



;Ui{X,t2) = iXui{X,t2). 



For the output yx{t2) = 
equivalent to the system o: 



?/l(A,t2) 
?/2(A,t2) _ 

equations 



, we shall obtain that fl2.17p is 



(A - 27rn(t2))?/i(A, ts) ~ iiB: + /?(t2))?/2(A, ^2) 



(/3(t2)-4)?/l(A,t2)+^2/2(A,t2 



from which we immediately obtain that ?/2(A, ^2) = i(/3(t2) — ^)l/i(A, ^2) 
and plugging it into the first equation 



^2/i(A,t2) - (7ril(t2) 



(3'{t2) + l3\t2))yi{X, t2) = tXyiiX, 



which means that yi{X,t2) satisfies the Sturm Liouville differential 
equation with the spectral parameter iX and the potential g(t2) = 
(vrii(t2) + /3'(t2) + /3'(t2)). 



The first equation fll.3l) considered for Hq 



^0" 



Hi' 



becomes 



-^7rii(t2) 







Hi' 



from where we conclude that 



(5.23) Hi' = -/3(t2), Hl^ - H^' = -tnn{t2). 

Let us consider the differential equation (15.41) . where we use the nota- 
Hi' Hi' 



tion Hi — j^2i j^22 
parameters, we shall obtain 



PHI' 



d TTll 

i-Hl' - PHI' + ^{Hl' 
£Hf + zvrni/o' 



. Substituting the expressions for the vessel 
= -Hl\ 



Hf) 



11 



0, 

Hl'-Hl\ 
Hi'. 



4-Hl' + r^iiHf + mf + ^Hf 

and consequently, using the formulas fl5.23p the second equation results 
in 

d 



(5.24) 



dt- 
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Together the first and the fourth equations give 



Hi 



12 



21 ^ 



Ai/22^^^^^^12^^^22^^^21 

at2 



from where we obtain using f l5.23p 



(5.25) 



11 



^ = C eC, 



22 



Additionally, i^o has to satisfy Hq = ^H^oi. Using this relation we 
shall obtain that 



rrll 7W^12 
rT21 tt22 

K = 
K = 





'01" 






1 







' 


1 " 




1 






from where we conclude that 



iH, 



22\* 
) 1 



12 




12\* 
) ' 



H, 



iH, 



21\* 
) ■ 



(5.26) 

As we can see H^^ = h^^ is a real valued (arbitrary at this stage) 
function, and Hn is as follows 



(5.27) 
where r G 



Hn 



r — ZTTu 



1^21 



r + zvTii 



Let us perform the same calculations for Hi. From the algebraic equa- 
tion fl53|l 

Hiai^ + ai'Hl = -H,ai'H; 

we obtain that 

Hl^ + [Hiy Hl^ + {HfY 
Hf + {Hiy Hf + {HfY 



-Ho<7i ^Hq 



From the equation fl5.4p with z = we obtain as before that 
(5.28) 

(5.29) 



Ht^ - Ht^ = —H^^ + ^n^lH, 



Hi 



12 



11 ' 

'Hi' - ^Hl^ + PH, 



d 



' 

11 



dU 
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and the same equation (15. 4p with i = 1 produces similarly to the pre- 
vious case 

d 



(5.30) 



H^^ - 13 Hi' + i{Hl' - Hi 



t22\ 



0, 



d 



(5,31) -JM 



11 



22\ 



~i^^,Hl' + P{Hi^-Hr). 



Plugging (I5:28|) and into flCTj) . we shall obtain that H^^ = h, 

have to satisfy the following differential equation of the first order: 

d 



21 



or after cancellations 



(5.32) 



d 



2z{H',')' = —,H^' - 2f3—H, 



rll 



dU 



Inserting here the expressions for H^ appearing in (I5.27P we shall ob- 
tain that the real part of the last equality can be derived from the 
equation (I5.24p : 



TTll 



0. 



The imaginary part gives the following equation 

d 



dt,' 



dt' 



-(tTh/?) + /^Vrii - P Till- TT^i - Tx'l^. 



_ 7 

Suppose (see |M2] ) that there exists a function r such that /3 = 



Then using fl5^ tth 



and inserting these equations into the 



(^) • 



formula for -^h^^ we obtain that 

dt2 T ^ T 

Let us write down the formulas for the elements of Hi. 

d 

dt2 

^dt2 - 



TJ\2 _ d tt21 TJU I iQUll 



HI 



11 



Hr = i{i-Hi^-PHi^) 



^ 'Hl^ + Hf) 



dt2 



21\l 



dt 



-innHl^ + p{Hl^ - Hf 
d 

dt2 



22\ 



Hi' - 2P-i-Hl' 



(E29D 

(E2H1) 
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The last equation f l5.32p ^ is obtained from fl5.32l) by substituting the 
index at Hq by the index 1: H^^ . 

Finally, we obtain that in the general case i^j+i is derived from 
using a system of similar equations. 



(5.33) 



TTl2 



TT22 



d ( TTll I rT22 \ 



from where we see that Markov moments are defined up to initial con- 
ditions for no = 2 elements. Notice that = 1 m this case. 
Let us also demonstrate theorem 15. 71 for the Sturm Liouville parameters 
from definition ESI We will take the simplest case n = 1 and as a result 
the transfer function is 



1 



\ + z 



C{t2)B{t2)ai 



oo ^ 
i=0 



-zYC{t2)B{t2)a^. 



We have already seen (in (15.241) ) that 7* is necessarily of the form 

/3 i 



for a real valued function /3(t2) on I. In this case, = which means 
that the first moment is a multiple of the zero moment: Hi = —zHq. 
Since the vessel parameters are constant the differential ring TZ = C 
is trivial. Using the formulas developed in |M2] we shall obtain that 
r = exp f j3 satisfies 



k)T = 0. 



for k = \/—iz, which may be rewritten as a polynomial differential 
equation for /3, after inserting the formula for r = exp J /3 and multi- 
plying by r-^ 

5.6. Non-Linear Shrodinger equation parameters. 

Definition 5.10. Non-Linear Shrodinger equation parameters are given 
by 



CTl = 

7*(x) 



1 
1 



(^2 



/3(x) 
-P*{x) 



1 








-1 



7(x) 
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Then the output and output compatibihty conditions take the form of 
the classical non linear Schrodinger equation with the spectral param- 
eter iX 

d 

-u{x, A) = {iXA + Q{x))u{x, A), 



dx 



where 



ai, A 



i 
-i 



-ia2, Q{x) = -7*(x). 



We may perform the same calculations for the moments as in the 
Sturm-Liouville case. Using formula fll.2p and denoting n-th moment 



22 



isfy the following relations 



(5.34) 



we shall find that the moments has to sat- 



d 



while the first moment Hq is found from the linkage condition (12.1 5p 



^0 



and the entries , Hq"^ are found using two last equations of (I5.34p : 
^^{Hl')=P{t,)H^,' = P{t2W{t,), 



6. Inverse problem at infinity for given vessel 

parameters 

In order to discuss the reconstruction of the transfer function from its 
moments, we will show first that for almost arbitrary vessel parame- 
ters one can solve the hq moment equations, appearing in Theorem 15. 51 
Once a symmetric, equal to identity at infinity, function is constructed, 
we will use a Krein space realization for it. Finally, by a counterexam- 
ple, we will show that not all vessels parameters admit Hilbert space 
(with positive X) realizations. 
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We have seen that there is a differential equation f l5.17p . which connects 
the moments, and means that the next moments is roughly speaking 
a second derivative of the previous one. If we start to estimate its 
growth, we will end with factorial growth for the moments, because 
consecutive applying of derivative produces many terms due to Leibniz 
rule. Instead, we will assume that there is a "big" function M{t2) which 
majorizes all the moments. Of course such an assumption narrows the 
class of function, to which the realization theorem 16.11 is applicable. 
On the other hand, the algebraic equations flS.Sp connects the real/ 
imaginary parts of the moments with all the previous ones. Again, we 
found assumptions, which will enable us to handle the growth of such 
elements. 

From theorem 15.21 it follows that we have to demand the symmetry 
condition (15. 5 p only at t^- Examining the equation (15. 5 p and supposing 
that ||i/n(t2)ll < fnC""^^, we will obtain that 

n 

lli7„W + (-i)vr'^:+iii <$^ii//n-,iiiiar^iiii//;ii < 

j=Q 

n n 
j=0 j=0 

So, if the series {/„} has the property 

n 

(6-1) fn+l = IWl^W'^fn-jfj, 

j=0 

we will obtain that 

n 
3=0 

One can find generating function F{x) = fix\ Multiplying (16. ip by 
x"^^ and summing, we will obtain that 



xF{xy = F{x) - F(0) F{x) = 1 4xF(0) ^ y. 

2x ^-^ 



oo 



n=0 



from where it follows that for -F(O) = 1 

4" 



nn{2n — 1) 
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Let M{t2) be a positive function with positive derivatives on [^2,^2], 
which satisfies the following differential equation 

, , C d , , 

= ||Ko||M(t2) + \\K^\\^M{t2) + WK^W^Mih) + ||K3||^M(t2) 
then, if Ef^{t2) < fnC'^^M^^^ti) we will obtain that 

II -TT -^n+l II ^ 

dt2 

</„C"+l(||iro||||i^n|| + ||i^l||||^i/n|| + ||/^2||||^/^n|| + ||^ 
</„C'"+l(||Ko||M(t2) + ||i^l||^M(t2) + ||i^2||^M(t2) + ||i^3||^M 

= -^C"+^M'(t2) < /„+lC"+2M'(t2) 

and by induction on further differentiating the equation fl5.17p and 

d^ 

using triangle inequality and additional assumption ||i^^j^^''(t2)|| < ^^||-^i(^2)|| 
we will obtain that 

Using these ideas we obtain the following 



^i5i(i2)|| < -^C"+2M(^)(t2) < /„,+iC"+2MW(t2) 



Theorem 6.1. Suppose that the vessels parameters and the norm are 
chosen so that Kq, Ki, K2, K3 satisfy for all k > 1 

(6.3) lli^f (t2)|| <-^||i^.(t2)||, z = 0,1,2,3. 

Suppose also that M{t2) is a solution of the differential equation (16. 2p 
with the initial value I, which majorizes i^o(^2) o,nd all its derivatives: 
II -^o'^^ (^2) II < foCM{t2) on [^25^2]- Then if at each step the imaginary 
(for even n) and real (for odd n) part of Hn+iCTi^ , which is uniquely 
determined from the first line in (15.161) also satisfies 

(6.4) \\H^^i{A)a^' - (-l)Vr^iJ:+i(t°)|| < 
then the series 

00 ^ 
^ - E A^^"(^2)ai 

n=0 
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determines an analytic at X = oo function S{X, ^2) iT'Ot necessary (11(^2)- 
contractive function (i.e. condition (12.191) of proposition [^751 does not 
necessary hold), which maps solutions of fl2.16p to solutions of fl2.17p 
and which is symmetric (i.e. satisfies fl2.23p ). Radius of convergence 

around X = 00 is at least ——. 

AC 

Proof: Let us write down all the ideas, preceding the claim of this 
theorem. We will prove the theorem, using induction on n. The basis 
of the induction follows from the assumptions, namely. 

Ili^ril </oCM(t2). 

Suppose that ^(ta)!! < fiC'^^ M^^\t2) for z = 0, 1, . . . , n, then from 
(15T7D it follows that 

II ^ ^ IK 

</nC"+l(||i^0||||i^,J| + ||/^l||||^i/n|| + ||i^2||||^//n|| + ||i^3|^ 

2 2 2 

</nC"+l(||i^0||M(t2) + ||i^l||^M(t2) + ||i^2||^M(t2) + ||X3||^M(t2)) 

2 2 2 



3^2 

and by induction on k, we obtain that for A; > 1 

II^Sl(^2)|| = 

k 



C"+^M'(t2) < /„+iC"+^M'(t2) 



using assumption 16.31 



, Z 

i=0 



' ^k 



i=0 



f^C-+\^M''f-'^ < fn+,C-+'M^'\t2). 



< 
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n 

3=0 

n n 

i=o j=o 
From (16.41) and this inequality it follows that 

\\H^+^{tl)a^'\\ < i||i/„+i(t^)ar^ + (-l)Vr^i/:+i(t^)|| + 

+ l\\H^+iif,)a^' - (-i)vr^i/:+i(t^)ii < 

Then 

Il^n+i(t2)ar1 = ||//„+i(t^)ar^ + r \\K_,M\\dy < 

Jtl 

Jq 

= /„+iC"+2 + /„+iC"+2(M(t2) - /) = /„+iC"+2M(t2). 
From these estimates we obtain that 

and clS cl result its radius of convergence is 

4 = lim ^/Hn{t)2 < lim y^^C^+^M^ = AC. 

XI n— >oo n— >oo 

because /„ ~ and the proof is finished. □ 



Using results in |DLdeS] we can easily show using Cayley transform 
that there exists a Krein space realization for the just constructed func- 
tion S{\, t2). 

Theorem 6.2. Let S{X, ^2) be an analytic at infinity function with 
value I there symmetric with respect to the imaginary axis (i.e. satisfies 
f l2.23p ), then there exists a realization 

S{X, t2) = I- B{t2rX-\t2){XI - A,)-'B{t2)ai{t2) 

such that the Lyapunov equation holds 

AX(t2) + Mt2)Al + B{t2yai{t2)B{t2) = 0. 
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Proof: The idea behind this theorem is that one can construct a Krein 
space reahzation at the given time and then obtain a vessel with an 
inner Krein space /C instead of "H, using the formulas in theorem 14.11 
Indeed, suppose that there is constructed a Krein space realization at 
t\ (which will be done in (16. 8p ): 

^(A, tl) = I- B;X,\\I - A^)-'Boa^i4). 

Define next B{t2) as the solution of (14. ip 

^[B{t2)a,{t2)] + A^B{t2)a2ih) + ^(^2)7(^2) = 0, B{4) = B^ 
0X2 

and X(f2) as a solution of (14. 3 p 

-^x(t2) = B{t2)<y2{t2)B{t2y, X(tO) = Xo. 
0X2 

and define '~i\{t2) using the Linkage condition (I2.15P for which we shall 
denote by $i(A,t2,i^2) the fundamental matrix of (12.171) . Then we will 
obtain that 

Y{\M) = ^l{\t2A)s{\A)^{\MA)- 

On the other hand, in theorem 16.11 we proved that one can construct a 
function S'(A,t2) with suitably chosen initial values, i.e. with S{\,t^ 
such that 

s{\t2) = ^A\t2A)s{\Am\t2A) 

for the given 7*(t2)- It is remained to notice that 

Y{\t2)S-\\t2) = $i(A,t2,t^)$*(A,t2,t^) 
is entire and identity at infinity, i.e. is /, which means that 

Y{XM) = S{\M) 

and from here immediately follows (as at the second part of 14. ip that 
l*it2) = jlih) and we have obtained a Krein space realization for the 
function S'(A,t2) constructed in theorem 16.11 

It remains to show that there always exists a realization at t^- Define 
a function Q{X) using Cayley transform, which satisfies 

(6.5) 5(-a, 4) = (/ + '-Q{\)a,A)){I - '-Q{X)a,A)r\ 

Actually, this function is given by 

Q(A) = 2ta^'A){I - 5(-zA, t°))(/ + S{-zX, 4))'' 
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and is well-defined at the neighborhood of infinity with value there. 
Then from a simple equality, resulting from the symmetry condition 
fl2.23p considered with —iX instead of A 

= -{! + S*{-tX, 4))cr^\tl){I - S{-zX, 

it follows that Q{X)* = Q{X) and Q{X,t2) is zero at the neighborhood 
of A = oo. Thus |DLdeSt Theorem 3] Q{X,t2) admits the following 
Krein space realization 

g(A) = r+(I-A/)-^r 

where for a Krein space /C, there is a bounded self-adjoint operator A G 
L(/C), and a bounded F G L{S,)C). So, we can insert this realization 
formula into (16. 5p and simplify: 

S{-zX,tl) ={I+ '-Q{X)a,{tl)){I - '-Q{X)a,{tl))-^ = 

= {2I-I+ '-Q{X)a,{tl)){I - '-Q{X)a,{tl))-^ = 

= -I + 2{I-'-Q{X)a,{4))-^ = 

= -I + 2(1 - '-T+{A- XI)-'Ta,itt))-' 

There is a simple formula |BGR] for evaluating the inverse of a matrix 
in a realized form: 

{I-'-T^A-XI)-'Ta,{4))-' = / + ^r+(Ix - XI)-'Ta^{4), 

where A^ = A — ■^ai{t^)T^ . So, the last formula becomes 

S{-iX, h) =-I + 2{I - zr+(I - A/)-ir(Ti(t°))-^ = 
(6 6) = -/ + 2(/ + ^r+(Ix - A/)-irai(t°)) = 

= / + ir+(Ax - A/)-irai(tO) = 
= / - r+(ax - 2A/)-irai(tO) 

Let us define Ai = —iA^ then we obtain that 

Ai+A+ = -ilx +2(Ix)+ = 

(6.7) = -^(A - ^rai(tO)r+) + + '-Ta,itl)T+) = 

since A is selfadjoint. In the Krein space, the Krein space adjoint may 
be represented by an invertible selfadjoint operator X in (the Hilbert 
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space sense), when one considers the Krein space with its inner product 
■ , ■ as a Hilbert space with its initial norm (■,■): 

[m, v] = {Xu, v) , Vm, V E )C. 
Then the formulas for Krein space adjoint becomes: 

At = x-M*x, r+ = xr* 

and the last equation l\6.7h 



A, + At 



-Tahiti: 



r+ ^ 



Ai + x-M*x = rai(tO)r*x ^ 
AiX-i + x-iy4* = rai(tO)r* 

which is the Lyapunov equation (\3.8\j at after defining Bq = T and 
Xq = X~^. Notice that Lyapunov equation (13.81) will also hold for all 
t2 as a result of lemma [27il Moreover, substituting A for —iX in fl6.6p . 
we also obtain 



(6i 



S{X,4 



I - B;X,\\I - A^)-'Boa^{4) 



□ 

It turns out that there are examples for which there is no a realiza- 
tion with strictly positive X. For, example considering the Non-Linear 
Schrodinger equation parameters, defined in section 15.61 for constant (3 



7* 



(5 
15* 



we can easily find that the fundamental matrices for the input and 
output LDEs are (for convenience t\ = 0) 



62 



t2 







where 



~2F 



E 



+ 



E- 



r + - 



2T 



[E 





E- 



[E-E- 



fr-f)i?^ (£_ 



-A , 



2T 



2T 



-^/3*, E = exp(rt2). 



Let us search for a constant in t2 matrix S{X, 0) 



a(A) 6(A) 
c(A) (i(A) 



such that the constructed S{X,t2) = $*(A, t2)'5'(A, 0)$ ^{X,t2) satisfies 
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conditions of the proposition 12.51 Considering the identity at infinity 
requirement, we come to the conclusion that 

c(A) = ^a{X), b{\) = -^d{X), 

where the square root at the Right Half Plane (RHP) is taken so that 
at the infinity its argument coincides with the argument of ^- this 

case the expression F + | is analytic in the whole C, except for a cut 
from —a/2 1/3 1 to a/2|/3|. In other words this can not be defined as an 
analytic function at the right half plane. On the other hand, if it were 
possible to realize this function with Ai, satisfying Lyapunov equation 
(13.81) with positive X, this would mean that the spectrum of Ai is at 
the left hand side, since ai = I and is a positive matrix: 

AiX + XA* = -B*B < 0. 

As a result S{X, 0) would have spectrum only at the left hand side of 
the complex plane, contradicting the previous argument. 
Finally, we present in this section a lemma, which says that one can 
always define moments so that the corresponding Pick matrices are 
positive, but as the example above shows they will not always define 
an analytic function (the radius of convergence may be zero) 

Lemma 6.3. Suppose that equations (15. 5p . (15. 4|) hold so that the eigen- 
values of the Pick matrix P„ are strictly positive on a finite interval I. 
Suppose that initial values are chosen for H2n+i, then there exists an 
initial values for H2n+2 such that all eigenvalues of Pn+i are positive 
on arbitrary finite interval around tg- 

Proof: Let us use Sylvester criterion for positive definiteness of a 
matrix. Denote the entries of the matrix (—1)'^ H2nO'i^ as follows 



and suppose that P„ > on a fixed interval X including tg- We will 
consider main minors of the matrix P„+i. First minors come from P„ 
and result are positive, the next one is 



hu hi2 

^21 ^22 



h 



Ip 
^2p 



Dn = [di 



hrA h., 



p2 



h 



PP 
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where is the first column of the matric C„. Its determinant, using 
a formula for block matrix is 



det 



det P„ det (/in + rfn - (C7^)*P-^C,^^ 



In order to obtain here a positive result, we have to demand that (since 
det P„ > and since /in + dn is real valued) 



3?/in > {ClTK'Cl 



Using the continuity and taking /iii(t2) big enough we obtain that the 
same inequality will hold on the interval X. Consider now the next 
main minor 

Pn Cn 

{C'J* hu + dn 

{ClY /i21+rf21 

Its determinant can be expressed as 

P. 



hl2 + di2 
/l22 + (^22 



det 



det 



{ClY hn + dn hu + d^ 

{C^y h21+d21 h22 + d22 
Pn Cl 

{Cir hu + du 



X 



p 


cl ' 


-1 


Cl 


_ {cir 


/in + dx\ _ 




hi2 + di2 



X det (/i22 + d22-[ [ClY /i21 + d2l ] 

from where we obtain that the necessary condition for its positiveness 



IS 



3ft/i22 > [ {ClY h2l + d2l ] 



p 




-1 


A^2 


_ {ClY 


/in + dn _ 




hi2 + di2 



22 



And so on for all the real parts of diagonal elements ha. It remains to 
remark that the conditions (15.51) . (15.41) do not involve real parts of ha 
and they can be determined from the positive definiteness of Pn+i- D 



7. Generalized Nevanlinna-Pick interpolation problem 

We address the two kinds of Nevanlinna-Pick interpolation problems 
11.41 and II. 5[ defined in the introduction. The first one is solved using 
the notion of positive pairs and for the second one we present a criterion 
in terms of the of the initial data. 
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7.1. Linear Fractional Transformations in terms of intertwin- 
ing positive pairs. Suppose that we are given a data of the NP in- 
terpolation problem 11.41 Following the notations of corollary 13.31 let us 
write 



, Ai = diag[-u)i, . . . , -u-;]. 



-Ut2ys{wn,t2y ut2y J 

Let us denote by capital Greek letters the following vessel parameters 



S2(t2) = 

r(t2) 



^l(t2) 

ai(t2) 

cr2{t2) 
0^2(^2) 

' l*{t2) 

7(^2) 



J, 



then simple calculations show that ^(^2) satisfies (14. ip 

-^[S(t2)Si(t2)] + AS(t2)S2(t2) + B{t2)T{t2) = 0, B{4) = B. 

dt2 

Suppose that X(t2) > is a solution of 

yllX(f2) + ^{t2)Al + S(t2)Si(t2)5*(t2) = 0, 
-^X(t2) = S(t2)S2(t2)S*(t2), 

dt2 

which is always possible if '^Wi 7^ for each i = 1, . . . ,n. Then the 
following collection 

Qj = {A, 5(t2),x(t2);Si(t2),S2(t2),r(t2),r,(t2);C2"; 

is a vessel for T^{t2) defined from the linkage condition fl2.15p 

r*(t2) = r(t2) + Si(t2)5(t2)*x-i(t2)s(t2)S2(t2) - 

-S2(t2)5(t2)*X-l(t2)5(t2)Sl(t2). 

Transfer function of the vessel QJ is 

e(A, t2) = hp - B{t2yX-\t2yXIn - A,)-'B{t2)^,{t2), 

which is in 5X(U* © U, U) for 

u = AS2(t2)-Si(t2):^ + r,(t2). 

dt2 
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e(A,t2) 



eii(A,t2) ei2(A,t2) 

e2l(A,t2) e22(A,t2) 



then if one defines 5*0 (A, ^2) so that 

(en + Soe2i)-^(ei2 + 5oe22) = 5, 

then the function ^(A, ^2) usually does not intertwine solutions of LDEs 
with spectral parameter A. 
Instead, we define 



so that 



1^1612 + W2Q22 = s, 



then the following lemma holds 



Lemma 7.1. The pair of functions W{X,t2) is in «SX(U,U*) and 

— — ^ — ^' ^ > on the domain of analyticity of W(X,t2) ■ 

A + /i 

Proof: Since 0(A,t2) is invertible for all A out of the spectrum of Ai, 

yx{t2] 



an element of U is of the form 6(A, ^2) 



Ux{t2 



where yx{t2), ux{t2) 



satisfy (I2.17p and (12.161) respectively. Then 

w{x,t2)e{x,t2 

= [ W,{X,t2) W2{X,t2) ] ^ 

= [ WiQn + 1^2621 mei2 + 14^2622 ] 



eii(A,t2) ei2(A,t2) 
e2i(A,t2) e22(A,t2) 



yx{t2) 

Ux{t2) 

yx{t2) 

uxih) 

yx{t2) 

Ux{t2) 



[ I S{XM) ] 



yx{t2) 

uxih) 



yx{t2) + S{X,t2)ux{t2) e U*, 



since yxih) and S{X,t2)ux{t2) are in U^.. 

From the formula WQ = [ ^ ] , it follows that 



w{x,t2) = [I s{x,t2) ]e~\x,t2). 
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A + /i 



considered on 



the domain of analyticity of W{-,t2) becomes 



A + /X 

e-i(A,t2)Si(t2)e-i*(Ax,t2; 



I 



A + 

Since B(A,t2) is a transfer function of a conservative vessel, its inverse 
is a transfer function too and satisfies 

0-i(A,t2)Si(t2)e(^,t2)-'*-Si(t2) 



A + yU 



> 



and consequently, 

l^(A,t2)Si(t2)W^*(/i,t2) 



A- 

> 

> 



> 



1 



X + fi 



[ I S{X,t2) ]Si(t2) 



x + fl'- J L 

SiX,t2)a^{t2)S*{^i,t2)-CXi{t2) 



^[I S{X,t2)] 





(Tl{t2) 



I 

S*{t^,t2 
I 

S*{f,,t2) 



> 



A + ^ 



by the properties of transfer functions for vessels. 
As a consequence of this theorem, we define 



□ 



Definition 7.2. A pair of functions 

W{X,t2)= [ Wi{X,t2) W2{X,t2) ] 

such that the matrix W{X,t2) has full rank is called positive if the 
conditions of lemma \7.1\ hold: 

W{X, t2) e 5X(U, U,), W{X, t2)JW{X, t2) > on C+ 

Previous arguments result in the following "partial solution" of the 
Nevanlinna-Pick interpolation problem 11.41 

Corollary 7.3. There exists a one-to-one correspondence between so- 
lution of the Nevanlinna-Pick interpolation problem I.4 and the positive 
pairs W{X,t2), which are analytic at the interpolation points. 

7.2. Nevanlinna-Pick interpolation problem 11.51 Using the pre- 
vious results, we present a criterion for solvability of the Problem II. 5[ 

Theorem 7.4. Given C^^^ -valued functions ai, a2, 7, an interval 1 
and N quadruples < 4, ^j, rjj > where ^2 £ I, Wj G C_|_, ^j, rjj G C^^^ 
j = 1, . . . A^, and assume that the corresponding matrices Xj > 0. Then 
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there exists a solution of the Nevanlinna-Pick prohlem M.^ i.e. there 
exists a function S G IZSX satisfying S{wi,t2)Ci = Vi ^/ o,nd only 
if there exists n e N, matrices A^,X^,G C^""^)''^""^) with > 0, 
Bl e CP^("-i), Fi,- e C"^" such that for A, Bi, X^ defined by 



(7.1) B, 

(7.2) Xi 

(7.3) Ai 



Vi - 6 



X^ 





Xi 



A' 



Xi 



X. 



it holds that 

(1) A, = V,,A,Vrr\ 

(2) §{\I - A,)-^Bia^\ti)^'\\,tlt^)d\ = V,,B,. 
and the matrix X{t2), 

-12 



X(t2 



X,- 



B,{y)a2{y)my)rdy 



is invertible on the interval I. 



Proof: For each t^ all the functions which satisfy S{wi,t2)^i = Vi 
of the form Te,(^^(A, t^)), provided X, = ^!^i^^2)4* - ^.^^1(^2)^.* 



Here 



e.(A) 



/ + 



Vi Vi (^i 
X,(A + < 
^ ^iVi(^i 

Xi(A + 



vKi(^i 



> 0. 



Xi(A + 



X(A + w* 



and Sn € 7?.«S. Given a minimal realization 



^o(A) 



I - iBi,nx},r\xi 



J^,)-'Bla,{f, 



of 5*0, we shall obtain from formulas f l3.15p . f l3.16p . f l3.17p the formulas 
(17.11) , (17. 2p , (17. 3p in the theorem. In view of Theorem 14.81 a necessary 
and sufficient condition to obtain the same function S{\,t2) for every 
i is that there exist invertible constant matrices Vij such that the op- 
erators Ai are similar. In other words there must exist V^j such that 



A, 



V^,A^^'. 



Moreover, the second part of theorem 14.81 tells that 
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additionally the equality 

^{XI - A,)-'Bia^\ti)<^-\X,ti,ti)dX = V^jB^ 
must hold. □ 
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